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Abstract 

We present here a relation of different types of Friedrichs mod- 
els and their use in the description and comprehension of resonance 
phenomena. We first discuss the basic Priedrichs model and obtain 
its resonance in the case that this is simple or doubly degenerated. 
Next, we discuss the model with N levels and show how the probabil- 
ity amplitude has an oscillatory behavior. Two generalizations of the 
Priedrichs model are suitable to introduce resonance behavior in quan- 
tum field theory. We also discuss a discrete version of the Priedrichs 
model and also a resonant interaction between two systems both with 
continuous spectrum. In an Appendix, we review the mathematics of 
rigged Hilbert spaces. 



1 Introduction. 

The Friedrichs model is a model aimed to describe the basic features of res- 
onance phenomena. The basic idea is considering resonances associated to 
a Hamiltonian pair {Hq, H}, where Hq is the Hamiltonian for the "non per- 
turbed" dynamics. Hq has a simple non-degenerate absolutely continuous 
spectrum that coincides with the positive semiaxis. In addition, Hq has at 



1 



least an eigenvalue imbedded in the continuous spectrum. The total or "per- 
turbed" Hamiltonian has the form H = Hq + XV, where V is a potential and 
A a coupling parameter that it is usually taken to be real (to preserve the self 
adjointness of H) and positive. The potential depends on a form factor func- 
tion f{uj), which determines the existence and properties of the resonance. 
The action of the potential is to transform the bound state into a resonance, 
characterized by a point in the complex plane, as shall be described below. 
This point depends analytically on the coupling parameter A. This is the 
basic description of the model as originally introduced by Friedrichs in 1948 

m- 

The first step to show that the Fridrichs model is an excellent device in 
order to understand the machinery of decay in Quantum Mechanics accessible 
to physicists was given by Horwitz and Marchand [59]. After that, there 
were given several generalizations of the original model for various purposes 
including a description of unstable theory of fields. 

In the present review, we intend to discuss most of the known versions of 
The Friedrichs model together their applications to model various situations 
in which quantum decay appears. 

The Fridrichs model was conceived as mathematically rigorous and ex- 
actly solvable so that it could well serve as a toy model for a precise de- 
scription of quantum decay. Also its possible generalizations are enormous 
in number and vast in applications. The present review is a first step to 
collect these generalizations. In order not to make this paper excessively 
long, we have selected some of these generalizations and not included a few 
ones. Our selection has been biased by our own work in the field. Examples 
of generalizations of the Fridrichs model that we have not included in our 
review are: 

i. ) The Fano- Anderson model [171 121 EO]- This is a Fridrichs model in 
which the unperturbed Hamiltonian has a bounded absolutely continuous 
spectrum. This model is useful in solid state physics to analyze instabilities 
due to the presence of resonances (see references in [70]). 

ii. ) The Cascade model |33- This is a model for unstable field theory. 

iii. ) A typical example of a generalization of the Fridrichs model that 
cannot be solved unless we make some rather severe approximations is given 
in a description of the boson-fermion interaction in nuclei developed by our 
group IMlllQ]. 

We also do not intend to discuss some special physical features concerning 
to decay, such as the Zeno effect that deserves a whole monograph both by 
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its extension and importance [721 |S3]. Neither the possible relation of the 
Friedrichs model with a model for quantum systems with diagonal singularity 

[nmiiiM]. 

We are mainly concern in the study of resonances in the Friedrichs model 
ant its various generalizations and we have obtained these resonances by 
means of the resolvent and not through the S matrix. Consequently, we did 
not attempt to obtain in our examples the 5* matrix or the M0ller opera- 
tors, which being relevant in a study of scattering are not necessary for our 
purpose^. 

The list of topics to be discussed in the present monograph goes as follows: 
In Section 2, we introduce a the basic Friedrichs model having a simple and 
as well as a double pole resonance. We define basic features in a language 
which can be accessible to both mathematicians and physicists. We devote 
2.1, to the description of the resonance that emerges along the construction 
and properties of its corresponding Gamow vector, i.e., the state vector that 
decays exponentially. We also include a situation in which the resonance 
comes from a double pole of the resolvent. In this case, the exponential 
decay of the resonance is multiplied by a linear function of time. 

The original Friedrichs model considers that the unperturbed Hamilto- 
nian Hq has one bound state only. What if we assume that Hq has more than 
one bound states? Then, the solution to the problem becomes more complex. 
In Section 3, we consider the case in which Hq has two bound states which 
become resonances due to the interaction with the potential. This is the two 
level Friedrichs model. 

The study of the level Fridrichs model is of particular interest. Here, 
the decay behavior, given by the survival probability, although modulated by 
an exponential, becomes oscillatory. This fact has been observed experimen- 
tally |1U3] . which incidentally shows the interest of these particular model in 
physics. 

In Section 5, we introduce two generalizations of the Friedrichs model to 
describe resonances in relativistic quantum field theory. In the first case, two 
scalar relativistic quantum fields ip and (p with respective masses M and m 
have an interaction of the type ipf"^. This problem can be solved on a sector 
in which a particle ip decays into two particles (p. The solution is provided by 

"'^A study of the Friedrichs model from the point of view of scattering theory is given in 
|44) . Among the big number of papers and textbooks on the quantum theory of scattering 
let me quote some absolutely essential books [78l [TJ [20] . A quite interesting review of one 
dimensional quantum scattering is given in |32) . 



3 



a generalization of the Friedrichs model, which shows resonances if M > 2m 
and can be exactly solved. 

In the second case, we study the interaction between two boson fields. 
One is a local field and the other a bilocal field with a continuous bounded 
mass spectrum. The interaction between these two fields is quadratic. This 
situation can be solved by formulating it in terms of a generalized Friedrichs 
model that can be exactly solved due to the fact that the interaction in 
quadratic. 

We conclude this section with a discussion of a Friedrichs model suitable 
for virtual transitions, which is possibly the simplest second quantization of 
the Friedrichs model. It studies the formation of a photon cloud around an 
atom. 

In Section 6, we present a miscellaneous selection of Friedrichs models. 
The former is a discrete version. Here, the non-interaction hamiltonian is a 
harmonic oscillators plus a bath of harmonic oscillators. Then, an interaction 
is switched on between the first oscillator and the bath. The second one has 
the originality of producing an interaction between two Hamiltonians with 
continuous spectrum. One of these Hamiltonians, corresponding to an inter- 
nal channel, has a bounded continuous spectrum doubly degenerated. The 
second one, corresponding to the external channel, has a continuous spec- 
trum coinciding with the positive semiaxis and is infinitely degenerated. An 
interaction intertwines both Hamiltonians. As a consequence, the continu- 
ous spectrum of the Hamiltonian of the internal channel is transformed in 
a branch cut for the resolvent of the total Hamiltonian. This branch cut is 
interpreted as a type of generalized resonance similar to the one obtained in 
Yukawa type interactions [79] . 

We finish the paper with two Appendices. In the former, we discuss 
the concept and main properties of Rigged Hilbert Spaces (RHS). RHS are 
essential in the definition and presentation of properties of Gamow vectors, or 
vector states for resonances. In a second Appendix, we show how to calculate 
the reduced resolvent, which gives the localization of the resonances, for the 
basic Friedrichs model. 

Although we are not going to study in depth the physical consequences of 
this model, a comment could be in order here. With respect to the relation 
between branch points and van Hove singularities, it has been reported [51] a 
behavior which is contrary to which is obtained in the basic Friedrichs model: 
for low dimensional systems the van Hove singularity leads to a non-analytic 
dependence on the decay rate associated with the resonance on the coupling 
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constant. 

Next, we start with the definition of resonance that we shall use in the 
sequel. 

1.1 Definition of resonance: comments. 

A resonance can be defined into several ways. It can be defined through 
the resolvent, or the behavior of the analytic continuation of the S-matrix. 
Both ways can be shown to be equivalent in specific models like the standard 
Friedrichs model. 

Any definition of resonance in nonrelativistic Quantum Mechanics relies 
in the existence of two dynamics, a free dynamics, represented by the unper- 
turbed or free Hamiltonian Hq and a perturbed dynamics given by the total 
Hamiltonian H = Hq + V, where is a potential responsible of the resonance 
behavior. Both Hamiltonians are densely defined on an infinite dimensional 
Hilbert space T-L of pure states of the system under study. 

After the properties of the resolvent of a self adjoint operator, we now 
that the resolvents 

are analytic on the complex variable z with a branch cut that coincides with 
the continuous spectrum of Hq and H, that we usually assume to be equal, 
in both cases, to the positive semiaxis R"*". 
Next, we make the following assumptions: 

i. ) There is a dense set V of vectors in Ti such that both RhoU^, z) and 
Rh{iP, z) admit an analytic continuation through the cut. 

ii. ) For some ip ^ V, the partial resolvent Rui'ip^z) has an isolated 
singular point (in general a pole) at a point zq of analyticity of Rho{iP, z). 

Then, we say that the Hamiltonian pair has a resonance at the point zq 

This is the definition that we shall use in the Friedrichs model. A second 
definition of resonance is also very much used and is suitable for the Friedrichs 
model: 

Let us consider that the Hamiltonian pair {Hq, H} satisfies the sufficient 
conditions such that the 5* operator exists. In this case, an incoming free 
state ip (evolving under the free dynamics Hq) will undergo the interaction 
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given by the potential V := H — Hq and will be released as a free state y?. 
The relation between this incoming and outgoing free states is given by the 
S-operator in the form: 

if = Sip 

In the energy representation, we can write ip{E) = S{E)iIj{E). Thus, the 
^-operator is a function of the energy E (and eventually of other variables). 
In general, S{E) is a complex analytic function with no other singularities 
than a branch cut coinciding with the positive semiaxis (and we assume that 
]R+ is the continuous spectrum of both Hq and H [791 123] . 

If the function of complex variable S{z) determined by S{E) has analytic 
continuation through the cut (and this can be shown to be the case under 
very general causality conditions [72]) and this continuation has poles outside 
the negative part of the real axis (then these poles come in complex conjugate 
pairs), then these poles are assumed to be resonance poles [23] . 

Thus, each resonance pole has a real part Eji and a nonzero imaginary 
part —ir/2. The energy Eji corresponds to the maximum of the bump of 
the cross section (whenever this equivalence applies) or the resonance energy 
and 7 is the width of this bump (related to the mean life of the resonance) . 

The relation between the latter definition of resonance and other defini- 
tions (bump in the cross section, sudden change in the phase shift, etc) is 
discussed in [25] . 

There are other definitions of resonances in quantum mechanics that ap- 
pear in the literature. Under general conditions, these definitions are equiv- 
alent either to the resolvent or the S matrix point of view. These are: 

i.) Resonances are the eigenvalues of the dilated Hamiltonian. A dilation 
(or dilatation) is a transformation U{6) on L^(M^) depending on a complex 
parameter 9, defined as follows: 

f/(e)^(x) =e3'/'V'(e'x). 

Clearly if 6 is real, U{6) is unitary. Dilations are used to define a large class of 
nonlocal potentials called dilation analytic potentials [SS]. If is a dilation 
analytic potential, then define: 

H{e) := u{e)[Ho + v]u-\e) = e-'^'Ho + v{e) . 

Then, resonances are the solutions of the eigenvalue equation: 
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The value of the eigenvalue zji of H{6) does not depend on the parameter 6 
and these eigenvalues are resonances according to the definition which makes 
use of the partial resolvent. For details, see [88i|- Dilatations really provide 
a method to obtain resonances called the complex scaling method. 

ii. ) When the potential is spherically symmetric, one form of finding 
resonances is obtaining solutions of the Schrodinger equation with complex 
eigenvalues, xi^y ^o) satisfying the so called purely outgoing boundary condi- 
tions. This means that as k i — > oo, then x(r, zq) ~ e*'^'', i.e., only outgoing 
wave function can exists. These complex eigenvalues zq can be identified 
with resonances and are poles of the 5* matrix. 

iii. ) Resonances are also defined as complex eigenvalues of a dissipative 
operator in the context of the Lax-Phillips formalism. This can be found in 

There are relations between these definitions of resonances. As was men- 
tioned above, some of them are consequences of others. For a review in this 
subject, see [3j and references thereof. 



2 The standard Friedrichs model. 

We shall present in here the Friedrichs model in its most simple form, which 
has been described in many places [50l|59l[l5l|ia[38l|71[93[96l|85l[85 We 
can introduced it in two different languages: Hilbert space language or bra 
ket formalism. The latter makes use of the rigged Hilbert space construction 
that we describe in the last part of the present review. Although we shall use 
of the bra-ket formalism along our presentation, we discuss briefly the Hilbert 
space presentation first and then we shall go to the second one which is more 
familiar to most users. Then, in its simplest basic form, the Friedrichs model 
includes the following ingredients: 

1.- A (infinite dimensional) Hilbert space H. If we use the energy repre- 
sentation, this Hilbert space is given by: 

H:=C©L2(M+). (2.1) 

Here, C is the field of complex numbers, R"*" is the positive semiaxis M"*" = 
[0, oo) and //^(R"*") is the Hilbert space of square integrable functions on M"*". 
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We say that we are using the energy representation because the Hamiltoni- 
ans we shall use in this Friedrichs model will have a continuous spectrum 
coinciding with R"*", then the space of the wave functions in the energy rep- 
resentation is then L^(M+). 

After (12. ip . any ip ^l-i can be written in matricial form as: 

^ = ( " ) , (2.2) 

where a is an arbitrary complex number and ^{oj) an arbitrary function in 
L^(R+). The scalar product of two vectors in % is given by 




a*P+ / ^*{Lo)ri{Lo)dLO. (2.3) 



Here, /3 is a complex number and rjiuo) G 

2.- In order to create resonance phenomena, we need a pair of Hamilto- 
nians: a free Hamiltonian Hq and a total or perturbed Hamiltonian H, such 
that the pair {ifo? H} is able to create resonance^. We here construct such 
a pair. Let us begin with Hq. Its domain (the subspace of all vectors ip E H 
such that HqiP G Ti) is given by the vectors in the form (12. 2p such thallf] 
u(p{u) G L^(M+). The action of Hq on a vector in its domain is given by 

/ a \ ( Wo" \ 
^0 = , (2.4) 

V ^{yj) ) \ ) 

where is a real positive number, i.e., ujq > 0. Observe that the vector 

(2.5) 

is an eigenvector of Hq with eigenvalue cjq, i.e., -f^o|l) = i^o|l). Note that 



^See |88| for examples of pairs of resonant Hamiltonians. 

^by definition, the energy representation is those in which the Hamiltonian is diagonal, 
or for the case of the continuous spectrum without degeneracy, the multiplication operator. 
This is a consequence of the spectral theorem [55] . 
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) 




) 



(2.6) 



which means that Hq is the multiphcation operator, for the second compo- 
nent, on the interval = [0, oo). We conclude that, by construction, Hq 
has a non degenerate absolutely continuous spectrum that coincides with 
[0, oo), plus an eigenvalue, Uq. Since ujq > 0, this unique eigenvalue of Hq is 
embedded in the continuous spectrum of Hq. 

3.- The total Hamiltonian has the form H = Hq + XV, where A is a 
coupling constant and is a potential to be defined. In order to create 
resonances, V should produce an interaction between discrete and continuous 
parts of Hq. It is proposed in the form: 



where ip is the arbitrary vector of "H given by (12. ip and /(w) is a function 
on called the form factor. Observe that Vip G "H for each ip & H ii 
and only if f{u)) G L^(M+). If this were the case, the domain of is 
(and V is bounded!) and the domain of H coincides with the domain of Hq. 
This may suggest that we cannot allo\\0for f{oj) outside L^(R+). The total 
Hamiltonian H has a continuous nondegenerate spectrum equal to [0, oo). 

So far, we have introduced the basic features of the Friedrichs model, 
without a discussion of its properties. Let us translate the above definition 
in terms of the bra-ket formalism. 

We have already defined the ket |1) G For each u in the (absolutely) 
continuous spectrum of Hq, there is an (generalizecifl) eigenvector of Hq, \uj), 
with eigenvalue u, i.e., Ho\co) = co\u). Taking this into account, we can use 
the following notation for an arbitrary vector ip as in (12. 2p : 



'^However, the RHS formalism allows non square integrable form factors, the Hamilto- 
nian would be in such case a mapping from the space of test vectors $ into its antidual 
$ ^ . See last section for this terminology. 

'^It is generalized in the sense that it does not belong to the Hilbert space, but to a 
bigger space that can be defined in the context of rigged Hilbert spaces [H]. Nevertheless, 
physicists use these generalized eigenvectors in their calculations. 




(2.7) 




(2.8) 
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The expression for the free Hamiltonian can be derived from the spectral 
representation theorem [SB] (see the specific form to translate integral spec- 
tral representations into the bra-ket formalism in [10] and a generalization 
in [52]). In our case, it takes the following form: 

Ho = ujo\l){l\+ u\u){u\du, (2.9) 
Jo 

where |1)(1| is the projection into the bound state |1) of Hq and |cj)(w| a 
similar object that can be defined for \u'^. 

Apart from the coupling constant A, the interaction potential V can be 
written in the following form: 

roc 

V= [f{u)\uj){l\ + nco)\l){u\]du, (2.10) 
Jo 

where /(cu) is the form factor that determines the behavior of the interac- 
tion. We recall that the form factor should be kept square integrable on 
the positive semiaxis if we want to keep the model in the context of Hilbert 
space, although this restriction may be unnecessary. Also, the form factor 
could be real or complex, it does not matter. From (12. 9p . we clearly see that 
V produces an interaction between the continuous and discrete parts of Hq 
modelled by the form factor function. We want to repeat again that the 
interaction also depends on the form factor A, so that H = Hq + XV. 

We want to show that even this simple Friedrichs model shows resonance 
behavior. In fact, this resonance is produced when the action of the interac- 
tion "dissolves" the discrete spectrum of Hq in the continuum. This process 
depends crucially of the form factor f{uj). The position of the resonance will 
also depend on the coupling constant A [S]. 

In order to search for resonance behavior, we need to fix a definition 
of resonance and we shall use the definition proposed in [88] that we have 
presented earlier (see (11. ip ). Accordingly, let us consider the projection P : = 
|1)(1| into the subspace spanned by the eigenvector |1) of Hq. Pick ip G T-L, 
as in (12.81) . Then, Pip = a\l). Then, take the operator, called the reduced 
resolvent: 

P{H-zI)-^P, (2.11) 



^Observe that the spectral measure for Hq can be written as dE^^ = jw) (a;| dw. See 

m. 
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where / is the identity operator. Then, consider the function on the complex 
variable 2;, depending of the vector ip eH: 



R^{z) := {^\P{H - zI)-'P\^) = \a\' . (2.12) 



The real number \a\ is irrelevant in the discussion of the analytic properties 
of R^{z), which are only contained in the following function independent of 

(IItt^II)- (2-13) 



r]{z) H — zl 

If this function, has an analytic continuation and the continuation has 
an isolated pole at zq = Er — iV /2 with nonzero imaginary part, this pole 
should be associated to a resonance, according to the definition of resonances 
from the resolvent properties [88]. The real part of zq, Er is the resonance 
energy (or location of the maximum of the bump of the cross section) and F 
the width of the bump. 

By extension, fl2.13p is also called the reduced resolvent. 

Some sufficient conditions for this analytic continuation to exist and have 
the desired properties have been studied in [591 El- These are conditions 
on the form factor |/(u;)p. Under these conditions, we can find an explicit 
expression for rj{z) in fl2.3p 



T^{z) = z-uo-\^ r du . (2.14) 

Jo z — u 

Under the sufficient conditions mentioned earlier [H], r]{z) is an analytic 
function on the complex plane C with a cut coinciding with the positive 
semiaxis (the continuous spectrum of both Hq and H) and admits analytic 
continuations beyond the cut from above to below and from below to above. 
In the first case, the boundary values of the function on the positive semiaxis 
are given by: 

n+{x) =x-uo-X' r '^^"^^''.^ du (2.15) 
Jn X - u + iO 



and in the second: 



ri_{x) = X - coo - r doj . (2.16) 

X — CO — iO 
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These analytic continuations have respective zeroes at the points zq = 
Eji — iT/2 and its complex conjugate Zq = Er + iV /2 (observe that (12.1 5p and 
(I2.16P are complex conjugate of each other), under the mentioned sufficient 
conditions [H]. Correspondingly, the analytic continuation of the inverses 
ri^^{x) and ■r]Z^{x) have poles at zq and z^ respectively. 

It is important to remark that the meaning of the kets \uj), some of the 
above formulas like (12.8112. 10]) as well as the Gamow vectors that will be 
introduced from the next section onwards, have a meaning in the context of 
extensions of Hilbert space called rigged Hilbert spaces. We shall add in the 
appendices a discussion on this mathematical tools and their applications in 
the theory of resonances. 

2.1 Simple pole resonances. 

These resonance poles may be simple or multiple. The case of the Friedrichs 
model with simple pole resonances is the easiest to study, yet containing all 
the nontrivial information concerning resonance behavior. As we have seen 
before, these resonance poles come into complex conjugate pairs. 

The next objective is to find the Gamow vector for this resonance. With 
this goal, let us construct the general state vector in the following form: 

V' = a|l)+/ ip{co) \co) dco . (2.17) 
Jo 

We want to emphasize that discrete and continuous degrees of freedom 
are orthogonafl. This idea implies the following relations: 

(1|1) = 1, (l|w) = (wjl) = 0, {uj\uj') = 5{uj - u') . (2.18) 

The Gamow vectors are defined as the generalized eigenvectors of the to- 
tal Hamiltonian with complex eigenvalues. These eigenvalues coincide with 
the resonance poles (as poles of the analytic continuation of the resolvent). 
As resonance poles come into complex conjugate pairs (in the energy rep- 
resentation), let us denote by Zq and Zq the poles corresponding to a given 
resonance. Then, their corresponding eigenvectors (Gamow vectors) are de- 
noted by l/o), l/o), so that 

^Being given an observable A with both discrete and continuous spectrum, the space 
spanned by the bound states is always orthogonal to the space spanned by the scattering 
states. 
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H\fo) = zo\fo); H\fo) = z;\fo). (2.19) 

Thus, in order to obtain the Gamow vectors, we have to solve the above 
equations. The technique is the following: Let us write the eigenvalue equa- 
tion for all X > 0: 

=0. (2.20) 

Here \l/(x) is the (generalized) eigenvector of H with eigenvalue x. The set 
of eigenvectors \E'(x) can be written in the following form: 

roo 

*(x) = a{x) |1) + / (f{x, u) \u) du , (2.21) 
Jo 

for all X > 0. 

The next step is to carry f l2.2ip into f l2.20p . The result is: 



(if - x)^(x) = (wo - a;)a(x) |1) + A / f*{uj)(^{x,uj) du \l) 

Jo 

roo roo 

+ {co - x)(f{x,uj)\co) dco + Xa{x) f (co) \uj) dco = . (2.22) 
Jo Jo 

Taking into account that the bound state |1) and the continuum are 
linearly independent, since they are orthogonal degrees of freedom (see rela- 
tions f l2.18p ) and omitting the vectors |1) in the first equation and |ci;) in the 
second, relation fl2.22p becomes the following pair of equations: 



{coo - x) a{x) + X ^{x,uj) f*{u)duj = (2.23) 
Jo 

{u-x)<^{x,u) + Xa{x)f{u) = 0. (2.24) 
Equation (12.241) can be written as: 

ip{x,u) = X^^!^^^^^^ + c6{x-u), (2.25) 
X — CO ±iO 

where c is an arbitrary constant that can be chosen equal to one (c = 1), 
without any change in the essential results. We see that two solutions appear. 
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one for each continuation (from above to below and from below to above). 
Then, equation f l2.25p becomes 



{uo-x)a{x) + X^a{x) I ^^^^^.^ du + Xf*{x) = 0. (2.26) 







X — u ±iO 



If we denote as 



r/±(x) ■.= uo-x + X' I ^^^^ dco , (2.27) 

X — u ±iO 



equation fl2.26p yields 



a(x) = -X^^. (2.28) 



If we carry this results into fl2.22p . we finally obtain: 



?7±(x) 



|1) + A/ — 

X — uj ±tO 



(2.29) 



Let us take one of these two vectors, say \I'4-(x). This vector is a functional 
on a certain space of test vectors, as we shall discuss in the Appendix. 
When we apply '^+{x) to a test vector ip+ G the result gives a function 
on the variable x, (\l/+(x) which admits an analytic continuation on the 
lower half plane, (\E'+(z)|v?+). This continuation has a pole at the point zq, 
which is the zero of the function ri+{z). We assume that this pole is simple. 
Sufficient conditions for this kind of situations exists and are discussed in [ii] . 
Then, if we omit the arbitrary cp^ G we can write on a neighborhood of 

^+{z) = -^ + o{z), (2.30) 

Z- Zq 

where o{z) denotes a Taylor series on zq, which converges on a neighborhood 
of Zo (both C and o{z) are functionals on $+). 

Let us go back to (12. 6p . As the analytic continuation ^+(-2) exists, the 
eigenvalue formula fl2.6p should be also valid for any z at which '^+{z) makes 
sense, so that we can write 
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{H-z)<iJ+{z) = 0. 
Now, if we carry f l2.30p into ( I2.3ip . we get: 



(2.31) 



= {H- z)^+(z) = {H -z)C + {H-z) o{z) , (2.32) 

Z — Zq 

which gives 

[H ~ zo)C = HC = zoC . (2.33) 

This shows that the residue C in f l2.30p is precisely the Gamow vector cor- 
responding to the resonance pole Zq. In order to obtain the explicit form for 
this Gamow vector, we just proceed as one does in a standard calculus of 
residua. As the pole comes from a zero of ri^{z) ~ constant (2; — ^o) + o{z), 
on a neighborhood of zo, we have: 



constant 



Since [55l 



|1) + A 



1 



z — CO + iO 



+ o{z) . (2.34) 



Z- Zq 



z — uj + zq — uj + [zq — uj + iO)^ 



+ o{z) 



we have that 



(2.35) 



constant 



|1) + A 



zo~ u + 



\uj) duo 



+ o{z) . (2.36) 



Therefore, save for an irrelevant constant, we conclude that 



C=|/o) 



ID./ 

^0 



A/(u;) 



\uj) du . 



(2.37) 



Zo — u + iO 

With the vector \l'_(a;) as in (I2.29p . we can proceed analogously, now 
taking into account that the function ri_{x) has a simple pole at Zq, the 
complex conjugate of Zo- By repeating the same arguments, we can obtain 
the solution |/o) of the eigenvalue equation H\fo) = ZqI/o) as: 
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l/o) = |l)+/ -^^^^\oo)du. (2.38) 



g — W — iO 

So far the discussion of the simplest Friedrichs model. The presentation 
of the model with a double pole resonance comes next. 

2.2 Double pole resonances. 

In order to show that there exists a Friedrichs model allowing for double pole 
resonances, we have to construct an explicit form factor that produces this 
effect. From all the above discussions, we need to show that a given form 
factor produces a double zero for ri±{z). Let us choose /(cu) such that 



with 



P{uj) - (o) - I3)(ui - n , (2.40) 

where /3 is a complex numbei|f| (with nonvanishing imaginary part!). With 
this choice, we obtain the following expression for the function r]{z), defined 
as in (127[3]|27H!1 : 



,W^.o-.-.A={^-^(^-^)}. (2.41, 

Next, we propose a change of variables which moves our formulas from the 
energy representation to the momentum representation: 

2=/ . ^ = ^2 ^2.42) 

and write ip{p) = Ti{p'^) = i]{z). Obviously, p is a complex variable. We have 
that 



*Note that (|2.39p behaves as w '"^^^^ and is therefore integrable on [0, oo), so that f{oj) 



is square integrable. 
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We shall show later that 77(2;) has a double zero if and only if ip{p) has 
a double zero. This is the first idea to be considered in order to determine 
under which conditions our model can have a double pole resonance. In 
addition, we have to fix certain constants. First of all, we have the complex 
constant b. This gives two undetermined real constants. But we can also 
play with the real constants Uq, A and with the real and imaginary parts of 
Po, the double zero of (p{p). Should {p{p) have a double zero at the point po, 
it might fulfill the following conditions: 

^(po)=0 ; <^'(po) = ; ¥^"(^0)7^0, (2.44) 

which give four equations, given by the real and imaginary components of 
the two identities in (12.441) . and a relation for the six parameters mentioned 
in the previous paragraph: uq, A, Realpo; Impo; Real 6 and Imb. Taking into 
account that we have two equations and six parameters, we have in princi- 
ple two free parameters that can be chosen arbitrarily among these six. It 
seems natural to choose the real constants uo and A as free, since they are 
usually data in the Friedrichs model. Then, after not complicated algebraic 
manipulations, we have 



CO, 



1/2 



2i 



TT 



A^ 



1/3 



(2.45) 



Po 



ttA^ 



2/3' 



— I 



TT 



A2 



If we carry fl2:45|) and (1236]) into fl2^ . we obtain 



Vip) 



{p-poY{p + plf 



1/3 



(2.46) 



(2.47) 



(P + -r) (p 

where r is some complex number different from p^, which is irrelevant in our 
discussion because it corresponds to a simple zero of the reduced resolvent. 
Equation fl2.47p shows that (p{p) has a double zero at the point po and another 
double zero at its complex conjugate Pq. 

Now, we have to show that the zeroes of ri{z) and (p{p) coincide, after the 
change of variables, as well as their multiplicity. To begin with, let us assume 
that ri{z) has a double zero at Zq. Then, there exists a function x(z) analytic 
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on a neighborhood of Zq, with xi^o) 7^ 0; cind such that on this neighborhood 
one has: 

The function if{p) = fjip^) and therefore, 

If A{p) := {p + PoYxip'^), we have that A(po) = 2poX(Po) 7^ and 

ip(p) = {p-po)^A{p) 

and ip{p) has a double pole at pq = ^/zq. Conversely, let us assume that ip{p) 
has a double pole at pq. Then, there exists an analytic function C,{p) on a 
neighborhood of po, with C,{po) 7^ 0, such that 

ip(p) = {p- pof^ip) = v{z) ■ 

Then, 77(2;) has a double pole at Zq = p^ if and only if the function (^l^^]-)2 has 
a nonvanishing limit if z goes to zq. Thus, 



which proves our claim. The conclusion is that ri{zo) has a double zero at zq 
if and only if ip{p) has a double zero at po with zq = Pq. 

Consequently, the analytic continuation through IR+ of f l2.10p has a double 
pole at the points Zq and Zq and therefore the version of the Friedrichs model 
here discussed shows a double pole resonance. 

We would like to obtain the Gamow vectors for this resonance. Once we 
have established the existence of a double pole resonance, on a neighborhood 
of zq^ the vector '^+{z) as in (12.291) has the following behavioiQ: 

^+(^) = j-^, + + o{z) . (2.48) 

[Z - ZqY Z- Zq 

As in the case of a simple pole resonance, we have on a neighborhood of zq 



^As '^+{z) represents a functional on a certain space of test vectors (to be specified 
later) this behavior corresponds to (0+|^'+(z)) for each € 
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constant 
(z - z,Y 

constant 
{z - z,Y 

\{z - Zq) 



ii)+A r 

Jo 

ii)+A r 

Jo 



z — uj + iQ 
2o - w + zO 



\uj) duj 



(2;o - + ^0)^ 



uj) duj\ . 



(2.49) 



Then, 



Ci = ii)+A r 

^0 



Co 



-A 



2;o — + ^0 



/o (2:0 - a; + i0)2 
Now, taking into account f l2.34p . we have: 



\uj) doj 
uj) duj . 



(2.50) 
(2.51) 



{H-z) 



H h o(z) 



0. 



(2.52) 



Xz-zoY Z-Zq 

If we multiply (I2.52p by {z — zqY and take the limit as zq \ — > 0, we obtain: 

{H - ;zo) Ci = ^ HCi = zqCi , (2.53) 

since o{z) is a regular function on z. If now, we multiply fl2.49p by z — zo, 
we get 



(H-z) 



Z- Zq 



+ C2 + {z- zo)o{z) 



Zq- z 
z- Zo 



Ci + {H- z)C2 + {z- Zo){H - z)o{z) = 



Taking the limit as z 1 — > Zq in (12.540 . one readily obtains 

HC2 = Z^C2 + Ci . 



(2.54) 



(2.55) 
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Summarizing, we obtain the following relations 



and HC2 = Z0C2 + C1. 



(2.56) 



The vectors Ci and C2 are linearly independent (otherwise C2 would be 
an eigenvector of H with eigenvalue zo in clear contradiction of f l2.55p . These 
vectors belong to the antidual in the rigged Hilbert space C "H C 
(see Appendix), thus, Ci and C2 are functionals. The restriction of H on the 
two dimensional subspace of with basis {Ci, C2} has the following form 



which shows a block diagonal form. This kind of structure has been discussed 
in several contexts in [331 EH ESI EH E5] • Conventionally, Ci and C2 are called 
1/0) and 1/1) respectively [5| 

As with the Friedrichs model with a simple pole resonance, we can proceed 
analogously with the vector \l/_(x) as in fl2.26p . now taking into account that 
the function ri_{x) has a double pole at Zq, the complex conjugate of zq. 
By repeating the same arguments, we can obtain the solution |/o) of the 
eigenvalue equation i/|/o) = -^ol/o) in fl2.38p and a second vector 
which is given by: 



The vectors |/o) and belong to the antidual of a rigged Hilbert 
space C "H C and are therefore functionals on the space to be 
defined in the appendix. As in the previous case, the total Hamiltonian H 
can be extended into the antidual and we can also prove the foUwoing 
relations: 



Therefore, on the subspace spanned by the vectors |/o) and the total 
Hamiltonian has the following form: 




(2.57) 




(2.58) 



^l/o) = ^ol/o) 



i^l/i) = ^o|/i) + |/o). 



(2.59) 




(2.60) 
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We may also study the time evolution of the Gamow vectors |/o), 
l/o) and l/i) (sometimes also called Jordan-Gamow vectors, see [30|). As |/o) 
and l/i) in one side and |/o) and on the other belong to different spaces, 
we must study their time evolution separately. A thoroughly discussion has 
been presented elsewhere [U [38] . The final result is given by 

e-''"\fo) = e-^*^«|/o) (2.61) 

e-*^|A) = e-^*^«|/i)-^t|/o) (2.62) 

and 

e-'*^|7o) = e-^'^o\fo) (2.63) 

e-'*^|/i) = e-''^o\h) + tt\fo). (2.64) 

As for the case of simple pole resonances, the above time evolutions can 
be either valid for all values of time or just for the values t > (in f l2.61|2.62p ) 
and t < (in (I2.63|2.64l) ). depending on the construction of the spaces $± 
[21 [381 [HIES]. 

3 Friedrichs model with two bound states. 

In the present case, we start with an unperturbed Hamiltonian of the follow- 
ing form: 



= a;i|l)(l| + a;2|2)(2| + / uj\uj){uj\du 

Jo 

2 

+ Y,[f:iu)\u){z\ + Mu)\z){u\]du. (3.1) 

i=l 

The most general vector in the space of pure states should have the fol- 
lowing form: 
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= y2^iiE)\i) + (j){E,uj)\uj)duj. (3.2) 
Let us carry f l3.ip into fl3.2p to obtain the eigenvalue equation. One gets: 



{u — E)(j){E, u) \u) duj 



2 «c 

[H -E)m = y^^,{E){uj,-E)\{)+ / 

poo r j'oo 



From f l3.3p . we obtain the following pair of equations: 



(3.3) 



{co, - E)^,{E) + ^{E,u)f,{u)dco = 0, % = \,1. 



(3.4) 



-E)0(E,u;) + J]/;(a;)^,(E) = 0. 



i=l 



(3.5) 



It is rather straightforward to obtain the solution to this system of equa- 
tions. If we write <^{E,uj) in terms of (pi{E), we get: 



(p{E,u) = cS{u - E) -J2 



1=1 



CO- E + iO 



(3.6) 



If we carry (13. 6p into (13. 4p . we obtain: 

fi{uj)fkiuj) 



[UJi 



i=l 



u- E + iQ 



ME) dco = -cfi{E) . (3.7) 



Then, ( 13. 7p give a two equation system with two undetermined functions 
'^i{E), i = 1,2. This system can be rewritten as follows: 
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io-E + iO 



■ duj 



/1M/2M 
u-E + m 



duj 



-cfiiE), (3.8) 



(wi - E)^2{E) - 



I/2MP 
io-E + iO 



■ duj 



V2{E) - 



du 



VliE) 



/o oj-E + iO 

= -ch{E). (3.9) 

Note that the determinant of the coefficients is given by: 



A 



00 



'^l /o u;~E+iO 
-Jo u:-E+iO ^2 -Jo a;-E+»0 



(3.10) 



Resonance poles are complex solutions on the variable energy of the equation 
A = 0. The solutions of the above system of equations (13. 8113. 9p are given by 



'PiiE) 
and 



A 



U)2 — E 



00-E + iO 



duj^ h{E) + h{E) j^' 



/iM/2(c^) 
u-E + iO 
(3.11) 



du 



ME) = 



Ui — E 



u-E + iO 



du}jUE) + h{E) 



f2{co)f^iu) 



u-E + iO 
(3.12) 

As we note earlier, the search for resonance poles is the search for solutions 
in the energy of the equation A = 0. At the first sight, this equation looks 
difficult to solve. Therefore, we may try for an approximation that gives us 
an idea on the behavior of the resonance poles. This approximation is the 
following: taking into account that for i = 1,2 



du 



Iq u-E + iO JO 



u-E 



du-i7r\ME)\' 



(3.13) 
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where PV stands for principal value. If we assume that this principal value 
is negligible, the value of the determinant A is left as: 



A K (o), - B - 'ijr|/i(£;)p)(a)2 -E- i7^\h(E)\- 



uj-E + iQ J \Jo u-E + iO 

^ {ioi - E){uj2 -E)- 27r(a;i - E)\h{E)\^ - ^{uj^ - E)\f2{E)\^ = .(3.14) 

Note that the last approximation in fl3.14p neglects the terms of fourth order 
in fi{E). In the last line of f l3.14p appears the approximate form of the 
equation A = 0. The complex zeroes of this equation gives the resonance 
poles, which depend on the values of the form factors fi{E). In the simplest 
case, fi{E) could be approximated by constants and equation fl3.14p gives two 
complex solutions. However, this approximation may not be compatible with 
the negligibility of the principal value. Then, equation fl3.14p is in general 
much more complicated than a equation of second order in the variable E 
and may give even an infinite number of complex solutions. 



4 Friedrichs model with N levels: oscillating 
decay. 

In this section, we shall discuss a generalization of the previous one: Now, 
we have N discrete levels imbedded in the continuous spectrum of the to- 
tal Hamiltonian. In the present case, the Hamiltonian is a straightforward 
generalization of (13. ip . Here, H = Hq + XV with 

N 



J' poo 

Hq = 'y^uJk\k){k\ + / uj\uj) {uj\ doj 

Jo 



k=l 
N 



POO 

V = T fk{u;)[\k){u\ + \u){k\]du, (4.1) 

k=l "^0 

where \k) represents the state with energy Uk > 0, k = 1,2, . . . , N. We 
shall assume no degeneracy, i.e., Uk 7^ cok' H k k'. The form factors /fc(w) 
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are assumed to be square integrable, just as in the simplest case. Here, the 
Hilbert space of states is given hy n = ® L^{R+). Thus, relations (KTE^ 
read in the present case as 



{k\k') = 5kk', {u\k) = {k\u) = , {u\u') = 5{u - u') , (4.2) 

with k = 1,2, . . . , N. The situation we can expect is that, once the potential 
AV^ is switch on, the N levels become resonances [H]. We want to study 
the eigenvalue equation f l2.20p . in the present situation. The form of "^{x), 
solution of f l2.20p . is given by 

ip{x,u)\u) du . (4.3) 



PC 

k=l -^0 



Then, we insert (14. 3 p into (I2.20p to obtain a system of equations, as (I2.23p 
and ([221: 



{uk - x)'ijjk{x) + X fk{uj)i){x,uj) duj = (4.4) 
Jo 

N 

{u-x)^{x,u) + Xj2fki(^)Mx) = 0. (4.5) 

k=l 

We can eliminate the function 4j{x, u) in this system to get the following equa- 
tion for -0^(3^) (note that ?/'(x, cu) = —[J2k=ifk{^)'4'k{x)]/{<-^~x)~C6{u — x), 
C being an arbitrary constant): 

AT 

J2[G-%k'{x)Mx) = -CXfkix) , (4.6) 

k=l 

with 

[G-%k'ix) = (cok - x)5kk' - r fjMMA . (4.7) 

Jo ^~x 

In (14. 6p . C is an arbitrary constant. We have to provide conditions for 
the analyticity of the function Gkk'{x), when a; is a complex variable. The 
solution of (14. 6 p is 
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N 

M^) = -ex Yl ^fc^''^^ ± *0)/,,(x) . (4.8) 

k'=l 

When the specific conditions for the functions [G^^]kk'{x) to have an analytic 
continuation are satisfied one can show that this analytic continuation 
includes a branch cut that coincides with the spectrum of H which is purely 
continuous in the positive semiaxis M"*" = [0, oo). They have analytic contin- 
uation from above to below and from below to above through the cut, that 
obviously cannot coincide with the values of the function of the complex 
variable x on the lower and upper half planes respectively!^. The respective 
boundary values on M"*" are denoted with the signs — and + respectively. 
This solution and (14.51) give us the form of the function tlj{x,u) as 



ip{x, (jj) = C 



X'T.k'=ih{^)[G-%Ax)U'{x) 



5{x 



X — a; ± iO 



(4.9) 



Once we have obtained the solutions of equations (14. 4114. 5p . we have two 
solutions of the eigenvalue problem {H — E)'^{x) = that we shall denote 
as \l/±(a:). Alternatively and following the typical terminology of scattering 
theory, these solutions are called "in" (for +) and "out" (for — ), i.e., \Ef_ 



[X 



^1 



and \E'_(x) = ^'out [12] • These solutions are: 



^±(x) = |x) + A V fi{x)GUx ± ^0) ( r 1^) _ . (4.10) 

Here, Hq\x) = x\x). The \x) in (I4.10p comes from the delta term in (14. 9p . 
We have chosen the value C = 1, which gives the following normalization 
condition: 

($±(x)|^±(x')) = 5(x-x'). (4.11) 

The solutions of the eigenvalue problem {H — £')\E'(i?) = 0, for E G [0, oo), 
are complete in the sense that 

^°This is the reason for defining these analytic continuations on the second sheet of a 
Riemann surface. Note that the analytic structure of the continuation is what it matters 
and not the Riemann surface that wc can ignore. 
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/■oo 

/= / |^±(x))(^±(x)|rfx (4.12) 
Jo 

and diagonalice the total Hamiltonian: 

POO 

H= x|^±(x))(^±(x)|(ix. (4.13) 
Jo 

Note that these resuhs are vahd for any value of including = 1 [T2] . 
From here, we can obtain a series of formal expressions whose validity will be 
checked in the next subsection. Taking into account that, by hypothesis, the 
bound states {\k)} and the kets in the continuum, {jw)}, form a complete 
set, we have that 

I = \^\k)(k\ + / \uj){uj\duj. (4.14) 

Assuming that identities in fl4.12p and (14.141) are the same, we immediately 
find the following formal identitieJ"1: 



\k)=l |^±(x))(^±(x)|A;)rfx ; \uj) = I |^±(x))(^±(x)|w) c/x . 
Jo Jo 

(4.15) 

The expressions {^±{x)\k) and {^±{x)\u) should be the complex conjugates 
of {k\^!±{x)) and {uj\^±{x)) that can be determined from (H?TU]) and (H^ . 
The result can be easily obtained to be 



N 

{k\m^{x)) = -\Y,fi{^)Gki{x±iQ), (4.16) 
1=1 

{.j\^^{x)) = - X) - T ^^MtM^l^IlM . (4.17) 

k,l=l ^ 

In order to obtain explicit expressions of \k) and |a;) in terms of |\E'±(x)), we 
have to carry fl4.16p and fl4.17p into the first and second equations in (I4.15p 
respectively. The results are 

^^As ^±{x) and \uj) are functionals, the precise meaning of expressions like (^±(a;)|fc) 
or {'i>±{x)\uj) should be clarified. 
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1^) = dxfi{x)Gki{x±iO)\^^{x)) (4.18) 

1=1 

Ic) = |vl/^(a;)) - A Mu) r dx ^f^^^)GKi{^) . (4.19) 

Note that G\^i{x + zO) = Gk,i{x — iO), where the star denotes complex conju- 
gation. These inverse relations will be used below. 



4.1 Survival probability. 

We know that the vectors |\E'-|-(x)) are eigenvectors of the total Hamiltonian 
H with eigenvalue given by u. Then, time evolution for |\E'+(x)) is given by 

e-''" = e"'*" . (4.20) 

This allows us to obtain the time evolution, with respect to the total evolution 
H, of the eigenvectors of the free Hamiltonian Hq, \k). Let us apply the 
time evolution operator e~**^ to (14.181) . with sign +. We readily obtain the 
following expression: 



\k)t:=e-''"\k) = -XY, / duje~''^ fi{uj)Gki{^-iQ)\m+{x)). (4.21) 

1=1 -^0 

This vectors \k)t must have a general similar to (14.31) . Also, inserting (14.101) 
into (I4.2ip . we obtain that this form should be 

N N 

\k)t = J2Mt)\l) + >^J2 dufiiu)gkiitu,t)\u), (4.22) 
1=1 1=1 -^0 

where |/), / = 0,1,...,A^ are the eigenvectors of the free Hamiltonian Hq. 
The explicit form of the functions Aki{t) and gki{uj,t) is obtained through 
a cumbersome but straightforward calculation, for which the main steps are 
explicitly given in [1^. The final results are for Aki{t) 
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Figure 1: Contour of integration C. 



(4.23) 



and C is the contour depicted in Figure 1. For gki{uj,t), we have 



2m 



1 



(4.24) 



where the contour C is again the same as in (I4.23p . There exists the following 
relation between Akiit) and gkiiu.t): 



as it can be easily derived from fl4.23p and fl4.24p . 

In order to obtain the survival probability for the bound states of Hq, we 
just calculate the survival amplitude for these states. Let |$) be the most 
general linear superposition of the eigenvectors of Hq, it has the following 
form: 




(4.25) 
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N 



\<^)=J2ak\k). (4.26) 

k=l 

Obviously, time evolution for fl4.26p must have the following form 

N 

l'^'W) = 5Z«fc|fc)t- (4.27) 

k=l 

We are now in the position of computing the survival amplitude A{t) and 
hence the survival probability |A(t)p. For the survival amplitude, we have 

TV N 

A{t) = ($|<l>(t)) = J2 = I] akal^Akkit) . (4.28) 

k,k'=l k,k'=l 

There is a form to write fl4.28p as a sum of the contributions of the resonance 
poles plus the background. We need to change the contour C into Ci as 
depicted in Figure 2. Note that for this contour transformation we are using 
the idea that the analytic continuation from above to below is produced in 
the second sheet of a Riemann surface. With these ideas in mind, we obtain 
the following expression for ^^^/(t): 

Akkit) = - V ri,, e-*^^* + ^ [ du e-'^' Guk'{uj) . (4.29) 

Here the sum in j runs out the poles of Gkk'i^) in the forth quadrant of the 
complex plane (see Figure .). The term r^f^, is the residue of Gkyiuj) at the 
pole Zj-. 



rik' = ^es Gkk'i(^] 



(4.30) 



Formula fl4.29p can be written in a more explicit form, provided that we are 
assuming that the form factors /fc(w) are analytically continuable into the 
lower half plane. After some calculations (given in [12] ), one obtains: 

N 

j k,l,k',V=l "^"^i 

(4.31) 
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X X 
X X 



Figure 2: Contour of integration Ci. By x we represent the location of the 
resonance poles. 
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where the Nj are certain normahzation constants |12] . 

The form of the terms Ai^ki{t) and hence the form of the survival amph- 
tude is particularly interesting. It has an oscillating term plus another term 
called the background and that is represented by the integral in fl4.29p or 
fl4.3ip . The background plays an essential role in the behavior of the decay 
phenomena at very short and very long values of time t. This behavior has 
been extensively studied 113 [TB SHI ESI EZ] • 

In the simplest Fridrichs model, N = 1 and there is a simple pole reso- 
nance only located at zq. Then, the form of A{t) is simply 

A{t) = -r e-**"° + — [ du e"*^* G{uj) . (4.32) 

In this case, the survival probability has a term that decays exponen- 

tially as e"*^™^" and that corresponds to the behavior of resonance phenom- 
ena for intermediate times, i. e., neither too short or two large 09]. However, 
if > 2, the survival probability has an oscillating behavior for which the 
maximus decay as t goes to infinite. 



4.2 Two level model. 

The simplest case of the level model is of course N = 2. Also for simplicity, 
we have chosen the coupling constant X in H = Hq + V equal to one, A = 1 
This situation has already several interesting features and for this reason and 
also for its simplicity it is an example of obligatory study. 

First of all, it is important to choose the form factors adequately. The fol- 
lowing choice was proposed in the context of a study of decay for elementary 
particles |68] : 

/fcM = ^^, A; = 1,2, (4.33) 

^ + Pk 

where pk are positive numbers. These form factors permit to calculate the 
matrix elements G^^, (u) from (14.61) as 

ttA^ 1 
Pk + Pk' iV^ + tpk){V^ + ^pk') 

with k,k' = 1,2. The branch square root is defined so that a positive real 
number had a positive square root. 
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From fl4.34p obtaining the 2x2 matrix G{uj) is just a straightforward ma- 
nipulation. In order to obtain the resonances as poles of G{u) the technique 
will be to obtain these resonances as zeroes of the analytic continuation of 
the inverse of the determinant of the 2x2 matrix G{u). If we make the 
simple change of variables given by y/u = ix, we need to solve the following 
equation in x: 



(detG'(a;))- 



{ui + x''){x + pif 



TT 



A2 



2pi 



{U2 + X''){X + - 



0. 



2p2 



TX 



Pi + P2 



(4.35) 



This is an algebraic equation of degree 8 with real coefficients. Therefore, 
all its roots are either real or appear in complex conjugate pairs. If certain 
conditions on tUj and pi are fulfilled [12], G{uj) has no singularities except for 
a brunch cut coinciding with the positive semiaxis. Then, all the resonance 
poles are poles of the analytic continuation, or in the language of Riemann 
surfaces, they are located in the second Riemann sheet. Real singularities 
correspond to virtual poles and lie in the negative semiaxis (second sheet). 

In the case of weak coupling, i. e., when the coupling constant A is very 
small, two pairs of complex conjugate poles can be evaluated perturbatively. 
The result is given by: 



TtA^ [^,-ip,r TT^A^ 
Zj = Ui + 7 ^-T^ h 



^^2p, (u;, + p2)2 -(^ + ,p^.)^ 
1 



X 



{Uj-Uk){pi + p2)'^{y/uJJ + ipky ip'^j ^ + ipjY 
+ o(A6), j = 1,2, kj^j. (4.36) 

In this case, we can obtain approximate expressions for the real and imaginary 
parts of Zj and the result is given by 
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uj := Real Zj = w • + — — i- + o(A^) , j = 1, 2 ^ 



7,:=-Imz, = ^-^^ + o(A^), j = l,2. (4.37) 

l^i + Pj) 

Note that (14.351) is of eight degree, so that it must have eight solutions. We 
have already calculated four of them. The other four are real and correspond 
to virtual states [T2], also called anti-bound states. 

Then, ( 14.37^ permits us to write an approximate expression for the sur- 
vival amplitud A{t). If we neglect the background term, (14.371) yields [I2] 



A{t) ^ 

fc,fc'=l,2 j=l,2 

^ ^ V2p,(p,-^v^)=^v^ 

27rReal (aia;) e~*^^* \ , , . 

(Pl + P2)(Pj - ^^/UJ]){pl - «v^)(Wj - UJl) J 

(4.38) 

It is important to remark that both in (I4.36P and (I4.38P the term 1/ {uk — oJi) 
appears. Then, these expressions cannot be used in the case of the existence 
of degeneracy levels for Hq. 

Note that, although for some initial conditions like Oi = 1, 02 = 0, the 
survival amplitude (I4.38P does not show oscillations, these oscillations must 
be present in the terms of the order A^. For the lowest order A^, the survival 
probability can be written in the following form 

\A{t)\^ = I |ai|2e-^^* + |a2pe-^^*e-2*"f , (4.39) 

with 

UJl—UJo , ^ 

V ■= — . (4.40) 
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From fl4.39p . we can see the dependence of the survival probabihty with 
respect to the initial conditions. For oi = or 02 = 0, the survival probability 
behaves exponentially without oscillations, at this order of approximation. 
If ai ~ 02 the survival probability shows considerable oscillations. 



4.3 The iV-level model. 

An analysis can be carried out for the N level model in the weak coupling 
regime. We can use arbitrary form factors fk{uj) in (14.71) . Then, we can get 



N N N 

(det Giu))-' = Hiuk ~u)~X^J2 n - ^) + ^(^') ' (4-41) 

k=l k=l m^k 



where 



Uu)= I du'^^^l^. (4.42) 
' u' — u — zO 



The zeroes of (I4.4ip give the position of the resonances: 



Zk = Uk- X'^hki^^k) + o{\^) = Uk- ijk, j = l,2,...,N , (4.43) 

where, again, Uk = Real^;^ and 7^ = Im^;^. 

Then, in the first nontrivial order, which is A^, we have for the real and 
imaginary part of the resonance pole located at Zk the following expression: 

u}k = ujk , 7fc = vrA^ fki^k) ■ (4.44) 

It is not a surprise that, at first order, real parts of Zk do not depend on A, 
while imaginary parts depend on A^. This is exactly what happens in the 
two level model studied in the previous subsection, see (14.36^ . At this order 
of approximation, we can also calculate the partial resolvent G jl2] 

Gkk'i^j) = (ujk-uj- A^4fc/(a;) )"^ 5kk' + o(A^) , (4.45) 

the residues 

rik, = -6kk'hj + o(A^) (4.46) 
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and the survival amplitude 



N 

A{t) = J] |afe|2 e"''"^* A2/|K)t . (4.47) 

k=l 

The conclusion that one arrives is that, although the survival probability 
decays, the decaying mode can be strongly oscillatory, depending on the 
values of the parameters a^, k = 1,2, ... ,N (initial conditions). One study 
of particular forms of this decay can be seen in [12] . 

5 Some types of Friedrichs-like models and 
their applications to Quantum Field The- 
ory. 

More sophisticated versions of the Friedrichs model have been used to intro- 
duced models for unstable relativistic quantum fields exactly solvable on a 
sector or under a well determined simplification or approximation. We want 
to discuss a few models that have been introduced elsewhere [9l [U |6T1 [5l |6] 
that are very useful to give us an idea on the construction and behavior of 
relativistic resonances, not making use of the S matrix formalism (for a study 
on relativistic resonances in the context of the S matrix see [2S1 EI])- 

Along this section, we are going to present four different situations, which 
are somehow related with unstable quantum field theory. In the first one, we 
consider two local fields with a cubic interaction. We show that, although 
the model is in principle not exactly solvable, it can be solved in a sector in 
which a particle of one of the spices considered decays into two of the other 
spice. In the second one, the interaction is between a local quantum field 
and a bilocal quantum field with a continuous and bounded mass spectrum. 
We choose a quadratic interaction between these two fields, which makes 
the exact solution tractable through a rather complicated generalization of 
the Friedrichs model. Due to the complexity of this example, we shall limit 
ourselves to give here a summary of it. Readers interested in a deeper un- 
derstanding of this example can go to the original paper [9]. Finally, we 
propose an example of "second quantized Friedrichs model", the Friedrichs 
model with virtual transitions, which is probably the simplest case of these 
second quantized models. 
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5.1 First model. 



Assume that we have two real scalar relativistic quantum fields ip{x) and 
ip{x), where x is the fourth component space time, with respective masses 
given by m and M, coupled with the simplest cubic interaction. The Hamil- 
tonian of the system is given by 

H^Hm + HM + V, (5.1) 

where 



Hm ^ J dx (02 + {Vipf + m^ip^) (5.2) 
Hm ^ J dx (V-' + {V^f + M'V') (5.3) 

V ^ xj dxi){x)^'^{x) . (5.4) 

The dot means first derivative with respect to time. By boldface letters we 
denote three-dimensional vectors. Four- dimensional vectors in Minkowski 
space are denoted by roman style letters. The products of two four vectors 
as well as the scalar products of two three vectors are denoted by a dot. 
We use the standard metric (-I-, — , — , — ) of Minkowski space. For example: 
k ■ X — koXo — k • X. 

In terms of the creation and annihilation operators, the above formulae 
can be written in the interaction picture as: 
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Hm = j rfqp(q)6t(q)6(q) (5.5) 

Hm = j rfku;(k)at(k)a(k) (5.6) 

V/(t) ~ '^k(iqic?q2 J d^x[a\\s.) exp{ik ■ x) + a(k) exp(—ik ■ x)] 
X [&^(qi) exp(igi ■ x) + 6(qi) exp(— igi ■ x)] 

X [6"^(q2) exp(ig2 ■ a;) + 6(q2) exp(-ig2 ■ x)] . (5.7) 
The quantum field operators in fl5.2ti5.4p are given by 

^(x,*) = /d£|.t(k)axp(*..) + .(k)exp(-*..)] 

^(x.*) = /dq|(,t(q)exp(.,.x)+Kq)exp(-.,.x)l (6.8) 

with the Lorentz invariant measure 



The creation and annihilation operators in (15. 5115. 7p satisfy the usual com- 
mutation relations: 



[a(k),at(k')] = {2nf2u{k)S{k-k') 

[6(q),fet(q)] = (27r)32p(q)5(q - q') . (5.10) 

The first goal is to obtain an appropriate asymptotic dynamics for this 
situation. Note that scattering theory applies if asymptotic completeness 
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holds true between the free and the interacting fields [22] • If this assump- 
tion is not valid, one tries to find another solvable evolution which satisfies 
the asymptotic condition and reestablish scattering theory as a comparison 
between the interaction field and the redefined asymptotic field. A typical 
example is the Faddeev-Kulish |16] removal of infrared divergencies. 

After integration over the three dimensional space, i.e., with respect to 
d^x, on the right hand side of fl5.7p . we obtain eight terms of products of cre- 
ation and annihilation operators with t-dependent exponents accomplished 
with three dimensional 6 functions of momentum conservation. According 
to the Riemann-Lebuesgue lemma [87], the asymptotic behavior of Vi(t) is 
defined by the behavior of these t-dependent exponents in the integration 
domain. For example, one of the terms on the right hand side of fl5.7l) is 



dkdqi dq2 a^(k) 6(qi) 6(q2) exp(i[u;(k)-p(qi)-p(q2)]t)(27r)^ (5(k+qi+q2) 

(5.11) 



The asymptotic behavior of the integral ( 15. lip as t i — > ±oo is determined 
by the term [uj(k) — p(qi) — p(q2)]. Notice that 

w(k) +p(qi) +p(q2) >M + 2m>0. (5.12) 

k+qi+q2=0 

Therefore, the integral (15. lip goes rapidly to zero due to the fast oscillations 
(indeed to the Riemann-Lebesgue lemma). Let us consider another term 



dhdqidq2 o^(k)6(qi)6(q2)e 



i[i^(k)-p(qi)-p{q2)]t 



(27r)3(5(k-qi-q2). (5.13) 



As for the case of the previous integral (15. lip , the asymptotic behavior of 
(I5.13P is defined by the quantity 



A(k,qi)=a;(k)-p(q0-p(q2) ■ (5.14) 

k-qi-q2=0 

The values of A(k, qi) depend on the masses of the fields and ip. If M < 
2m, then A(k, qi) is strictly negative: 



A(k, qi) < 2{m^ + ^y^^ - (M^ + k^^^ < , 



(5.15) 
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To show flS.lSp . we obtain the maximum of A(k, qi) on qi for each k. This 
maximum is at qi = k/2 (it is certainly a maximum) and the value of 
A(k, qi) at the maximum is 

2(m2 + ^)i/2_(M2 + k2)V2. (5.16) 

Note that if M < 2m, f l5.16p is indeed smaller than zero. Thus, A(k, qi) 
has constant sign and therefore the Riemann-Lebesgue lemma applies. As a 
consequence, the integral f l5.13p decreases fast as t i — > ±00. 

In contrast, let us assume that M > 2m. Then the maximum (I5.13p is 
bigger than zero. However, for fixed k the term 



A(k, qi) = (M2 + k2)V2 _ (^2 ^ _ (^2 ^ _ j^)2)i/2 (5^^^) 

is obviously smaller than zero for high values of |qi|. Let us consider the 
manifold for which 



A(k,qi) = 0. (5.18) 

The integral f l5.13p over this region does not vanish asymptotically and 
the same behavior can be observed on its complex conjugate. In other words, 
by relinquishing the stability condition M < 2m, we obtain an unstable field 
theory where the asymptotic condition for scattering theory fails. Therefore, 
following a standard procedure in field theory, we redefine the asymptotic 
evolution in the interaction picture as follows: 

Ks = Ho + Vasit), (5.19) 

where Vas includes the slowly decreasing part of Vj. There is an ambiguity in 
the definition of Vas and we shall take it as the sum of f l5.13p and its complex 
conjugate, so that Vas is Hermitian. This is the simplest choice for Vas- 



Vasit) = A / (ik(iqi(iq2 [a'^(k)6(qi)6(q2)e 



«A{k,qi)t 



+ 



+ a(k)6t(qi)6t(q2)e-A(k,qO*](2^)3 ^(k _ - q^) 



A / ci=^x[^(+)(x)<^(-)2(x) + ^^~\x)^^+^\x)] , (5.20) 
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where A was defined in fl5.14p and '^^^\ denote the positive and negative 
frequency parts of (f. Same for ip. 

• The above comments have an obvious physical interpretation: The 
decay process given by 

if) I — > ^ + ^ 

is only possible if the mass M of the ?/;-particle is bigger than the mass 
2m corresponding to two y9-particles. 

• Renormalization. As we shall see below the has produces a ultraviolet 
divergence in the equation for its eigenstates. In order to remove this 
divergence we add to has the appropriate counterterm 



H,,,, = - / d'x5M\^^\x)^^-\x) = / dk --— a\k)a{k) , (5.21) 

where the mass renormalization SM"^ is of order A^. The appearance of 
ultraviolet counterterms is due to our choice of the asymptotic inter- 
action. We can of course introduce some smooth cut off in Vas, but it 
will involve additional parameters to the asymptotic states and it will 
generally break the relativistic invariance of the system. 

In order to obtain the desired formula for the asymptotic Hamiltonian, 
we add has and Hc.t.'- 

Has{t) = has{t)+H,,t,. (5.22) 

In the Schrodinger picture, the time dependent exponents e^^^^'^'i-^)* in 
equation fl5.20p disappear. Now, we shall consider the eigenstates of Has = 
-f^as(O), which is a Hermitian operator in the Fock space corresponding to 
the creation and annihilation operators satisfying fl5.10p . 

From now on, we shall work in the Schrodinger picture, in order to make 
calculations easier. In the Schrodinger picture the total Hamiltonian H will 
be the following: 



41 



Has = Hm + Hm + Hc.t. + X j dkdqi dq2 K(k)6(qi)6(q2) 

+a(k)6t(qi)6t(q2)] i2nf 5(k - qi - q^) • (5.23) 

We shall show how the above problem of unstable quantum field can be 
solved exactly on a sector. Of course, as the interaction is not quadratic, no 
exact solution is know for all sectors. This sector corresponds to the decay 
of one ^/'-particle into two 99-particles. The Hilbert space of this system is 
given by 

H:=H^®H^®H^. (5.24) 
The eigenspaces of the operator 

N = N^ + 2N^ (5.25) 

are invariant subspaces of Has- The operators A^^ and are the number 
operators of ip and (^-particles. The asymptotic interaction does not affect the 
vacuum state and the one V'-particle state, which are the only nondegenerate 
eigenspaces of A^. The first degenerate case is the eigenspace of with 
eigenvalue 2, which is precisely "H as given in fl5.23p . In the subspace f l5.23p 
the term Vas in Has produces transitions of two y^-particles state to one ip- 
particle and vice versa. We should mention that not every two <y9-particle 
state will mix with one V'-particle state, but only with those in which the 
two (y9-particles are in the S-wave state. 

The most general linear combination of two y^-particles in the S-wave 
state with one ^/'-particle, can be written as: 

<|.(E,k) = /(E,k)|k) + J dE'fiE,E\k)\E\k), (5.26) 

where |-E",k) is the two particle S*— wave and |k) is the one ip particle 
state. The objective is to solve the eigenvalue equation 

{Has -E)<I>{E, k)=0, (5.27) 

where Has was given in fl5.23p . The method to solve this equation consists 
in replacing fl5.26p into fl5.27p . Thus, we can obtain the coefficients f{E,k.) 
and f{E,E','k). 
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We are interested in obtaining these coefficients in the unstable case char- 
acterized by the condition M > 2m, where M and m are the respective 
masses of the and if particles. The method is tricky as it requieres the 
use of "counterterms" and a particular choice of a normalization, but it is 
explained in detail in [8]. The final result is given by 

and 



where 

is the partial resolvent of the Hamiltonian Has- This partial resolvent has 
an explicit expression [6]. The function r is defined on the halfiine E > 
(4m2 + k2)V2 by 



t{E' -k') = j rfqidq25(k - qi - q2)5(E - p(qi) - ^(qs)) 

- 41^^- A)"^''(^-'*'"^ + ''^'"^'• 
(5.3l) 

where 6{x) is the usual Heaviside step function, dq^ and w(k) were defined 
in (15.91) and A is an arbitrary constant. 

In fl5.29p there is an ambiguity which arises due to the singular denomina- 
tor {E' — E). This ambiguity also exists in f l5.28p although it is there hidden 
in the partial resolvent fl5.30p . Consequently, we shall define two solutions 
of eigenvalue problem (15.271) that correspond to the incoming and the outgo- 
ing solution. The incoming and outgoing solutions are defined as boundary 
functions of analytic functions from above {E + iO) and from below [E — iO) 
the real axis respectively: 
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^ ' ^ M ' / T- 2a;(k)r/(E ± zO, k) 



X 



|A^.(2.)3£.£'^,-l_5|£'.k)-|k)l} 



(5.32) 

This formula explicitly demonstrates that the energy spectrum of the 
asymptotic states in the eigenspace of N, as defined in (15. 25 p . corresponding 
to the eigenvalue 2, is absolutely continuous over the interval [Eq, oo) where 
Eq = {Am?' + k^)^/^. The isolated eigenvalue associated to the one ^/'-particle 
state |k) has been dissolved to the continuum due to the interaction. This 
reproduces somehow the situation that appears in the standard Friedrichs 
model discussed in Section 2. 

5.1.1 The Gamow vectors. 

In [6], it has been shown that the inverse partial resolvent ri{E,\<i) admits a 
continuation on the variable E, ?7+(i?,k), from above to below through the 
branch cut [Eq, oo) and has a complex zero zr on the lower halfplane. In 
order to find resonances, we have to solve the equation ri{E, k) = 0. The 
explicit form of //(-E, k) is 

r?(E,k)=u;(k)-E-A^(27r)3(E-a) / dE' ^ (5.33) 

Jeo [E - E){E - a) 

where t{E''^ — k^) was given in (I5.3ip . 

Note that the interaction Vas in (I5.20p depends on a multiplicative cou- 
pling constant A. When A = 0, the interaction is switched off, exactly as it 
happens in the ordinary Friedrichs model. It has been shown in [6j that in 
absence of interaction (A = 0), E = a;(k) and that for small values of A one 
has that the solution of the equation r]{E, k) = corresponding to the lower 
half plane should have the following form 

zr = w(k) + X^z, + O(A^) , (5.34) 
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where 



^1 = -- ( 1 - M^J • (5-35) 



57r 



Note that fl5.9p shows that = a;^(k) — k^, where M is the mass of the 
particle. Then, (15.341) and (I5.35P gives 

ZR = u{k) - ^ (k^ + zMVy/' , (5.36) 



where we have denoted 



Mr = — (1-^) . (5.37) 



in 



M2 



This zero has a complex conjugate on the analytic continuation from 
above to below. Equation (I5.36P demonstrates that the zero of ri{E, k) cor- 
responds to a particle with complex square mass 

= - iMT , (5.38) 

which is the remnant of a ■^/'-particle. The eigenstate f l5.32p as a function of 
E has a pole at the point E = and its residue in this pole is the eigenstate 
with complex energy. Up to an irrelevant normalization constant this state 
is: 

POO 1 

$^(^^,k) = Av^-(27r)M dE'- \E',k)-\k). (5.39) 

Jeo ^ - Zr 

This Gamow vector has a clear meaning as an antilinear functional on a 
suitably chosen space of test functions The test function space should be 
dense on the Hilbert space (I5.24p of the states of two v?-particles and one ip- 
particle, which is contained in the Fock space of the states of particles and 
^/'-particles. Following the spirit of the construction for simple non relativistic 
resonances, a simple choice for the test function space, dense in the Hilbert 
space (I5.24p . is given by 

$ := [{ul r}S)® Sim?)] © s{m?) . (5.40) 

We recall that 
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• "H^ is the space of Hardy functions on the lower half plane. 

• S" is the space of all real valued complex functions that are differentiable 
at all orders such that they as well as their derivatives vanish at the 
infinity faster than the inverse of any polynomial (the one dimensional 
Schwartz space). 

• S'(M^) represents the space of all complex functions on the three dimen- 
sional real space M.^ having the same properties as the Schwartz space 
S. 

Then, the rigged Hilbert space is a triplet of the form $ C H C 
where is the dual space of Then, Gamow vectors belong to this dual 

A typical function in $ has the form F[E, k) + g(k') where F{E, k) e 
{Hi nS)® S{M.^) and g{k') e ^(R^). The variables E and k in F{E,k) 
represents respectively the total energy and the sum of the momenta for the 
two y9-particles. The variable k' is the momentum of the V^-particle. 

The action of the functional ( KM on the function K{E, k", k') := F{E, k") 
+g{k.') is given by 



(ir(E,k",k')|$^(^ii,k)) = XV2.{2n)' / dE'- {F{EX)\E',k) 

JEo ^ - Zr 

-{g{k')\k) = \V2.i27rf TdE'^^ + gik). (5.41) 

JEo ^ - Zr 

For the fixed variable k, the function F[E, k) is a Hardy function and, there- 
fore, the integral term of f l5.4ip converges. In the usual Fock space notation, 
a typical element of can be written as 

F{EX)\EX)+9{'^W)- (5.42) 
The functional ^^{zr, k) has the following property: 

Has ^""{ZR, k) = ZR <l>^(zfi, k) . (5.43) 

It may be interesting to discuss the proof of this latter statement. We 
have seen that <l>(£',k) has a meromorphic extension from above to below 
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on the variable E with a pole at the point z/j. Lets us call this extension 
^(^(-Sjk). The residue of $(7(z,k) at zr gives the Gamow vector ^'^{zii,k). 
On a neighborhood of zr one has the following expansion in terms of the 
complex energy z: 

^c{z,k) = — - — ^^{zR,k) + regular terms. (5.44) 

z- Zr 

The meromorphic extension of $(-E',k) allows us to extend f l5.27p . so that, 
we have: 

Has^c{z,k) = ZR<^c{z,^), (5.45) 

for values of the complex energy z in the lower half plane. If we bring f l5.45p 
into fl5.44p . we get: 



Has^^{zR,k) , , ^<^^{zR,k) 

h regular terms = [z — zr + zr) h regular terms . 

Z — Zr Z — Zr 

(5.46) 

If we identify the pole terms in the left and right hand sides of fl5.46p . we 
finally arrive to (]5.42p . This concludes the presentation of the model. 

5.2 Resonances in the interaction between a local and 
a bilocal boson field. 

In this example, developed in [9], we shall consider an exactly solvable model 
for relativistic unstable quantum fields. This model considers a quadratic 
interaction between a local scalar field ip with a mass M and a bilocal scalar 
field i/j with a continuous mass spectrum (we shall explain this terminology 
soon). The interaction makes the local scalar field unstable. The interac- 
tion is exactly solvable and it can be worked out as a generalization of the 
Friedrichs model. 

It is very convenient to write the local scalar field y:>{x) with mass M in 
terms of the creation and annihilation operators as follows: 

¥p(x, t) = j dk [a){k) e*^-" + a(k) e'^'^'" ] . (5.47) 

Note that k = (fco,k) and x = (a:;o,x). By boldface letters we always de- 
note three dimensional vectors. Four dimensional vectors in in Minkowski 
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space are denoted by Roman style letters. The products of two vectors in 
Minkowski space are characterized by a dot. In Minkowski space we use the 

metric H . For example, k ■ x = koXQ — kx. The Lorentz invariant 

measure dk. was defined in (15. 9p . 

The creation and annihilation operators in (I5.47p satisfy the following 
commutation relations: 

[a(k),at(k)] = (27r)3 2u;((k)(5(k-k'). (5.48) 
The Hamiltonian for the field ^{x) is given by 

= j dkuiik) a\k) a{k) . (5.49) 

Therefore, the field ip{x) is a free Klein-Gordon field. 

Next, we consider a bilocal scalar field i'ix^, q) with continuous mass 
spectrum and we shall describe this field as follows. Apart from the depen- 
dence of ip on the four vector x with components x^, = 0, 1, 2, 3, it depends 
on an additional real variable q representing an internal degree of freedom. 

For simplicity in our calculations and because it does not interfere in our 
discussion, we shall assume that ipix^j,, q) is an even function of q. This means 
that 



i(j{Xf„q) = ij{x^,-q) . (5.50) 

The mass operator affects to this internal variable and in its simplest form 
has the form 

= 4m - — . (5.51) 
oq^ 

Then, the field ip{x^, q) satisfies the following generalized Klein-Gordon equa- 
tion: 

{D-M^)i;{x^,q) = 0, (5.52) 

where □ is the usual D'Alembert operator. 

Bilocal fields were introduced by Yukawa |102] and Markov [71] in con- 
nection with extended particles. The motivation of its use here is described 
in [6] in which the solution ip{x^, q) of (I5.52p can be seen as the state of 
two particles. The fact that the second field is bilocal with a mass operator 
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with continuous spectrum is causes a quadratic interaction between these 
two fields to have non trivial features such as instabilities. 

The solution of equation f l5.52p with condition f l5.50p can be written in 
the following form 



where 

K={E^-k^- 4m2)i/2 = ^(^^) . (5.54) 

We change the variables in fl5.53p . in order to use as a new independent 
variable instead of k: 

dn dE 

- = (5,55) 
After (15.551) . equation (I5.53P reads 



i^{x,,q) = j^E j '^^^^^{B\\.,E)e-'+B%E)e—'). (5.56) 

Following the routine of second quantization, we replace the function 
i?(k, E) by an operator, that we also will call -B(k, E) for simplicity, and its 
complex conjugate -B*(k, E) by the adjoint operator -B^(k, E), satisfying the 
following commutation relations: 

[i?(k, E), i?t(k, E)] = (27r)^ K{k^) 6\k, - k'^) . (5.57) 
Then, the function ip{x^,q) becomes the operator 



Hx„q) = jyE J ^-^^^^{B\k,E)e-UB{KE)e—'). (5.58) 

Then, we have constructed the quantum bilocal field. Recall that the solution 
ilj{Xfj_,q) to our bilocal field, represents somehow the state of two particles. 

Now, let us assume that the fields ( I5.58P and f l5.47p interact and that the 
interaction is given by 
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Hint = -X / (i^x / dqilj{Xf,,q) f{q)ip{x^) . (5.59) 

J J —oo 

The function /(g) is called the form factor and we choose it to be a smooth 
even function. As always, the coupling constant A is assumed real and posi- 
tive. As we mentioned already, the bilocal field results as an approximation 
of the behavior of two relativistic interacting particles. The interaction of 
the bilocal field with the field <f{x) represents the interaction of the two for- 
mer particles with a third of fixed mass M. This is therefore, a model that 
approaches the behavior of three interacting particles, just like the other one 
described in the precedent section. 

If a{y) is the Fourier transform of f{q) and a(k), a^Ck) are the respec- 
tive annihilation and creation operators for the local field ip{x), the total 
Hamiltonian is given by 



(5.60) 

and the three momentum is: 



(5.61) 



where w(k) was given in ( 15. 9p . 

The next step is to diagonalize the four momentum f l5.60p and f l5.6ip . This 
means that we are looking for creation, 6^(ii^,k), and annihilation, 6(£',k), 
operators such that the four momentum components can be written as: 

r d^kdE 
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with n{E, k) = K{k^) as given in fl5.54p . To get this diagonahzation, we pose 
the eigenvalue equation 

[P^,b\E,k)]=k,b\E,k), (5.63) 

where = (-E, k) and the are given by fl5.60p . fl5.6ip . To solve (I5.63p . 
i.e., to obtain the creation operators b^{E, k), we make the following Ansatz: 

b^{E,k) = j dE' {T{E,E',k)B\E',k) + R{E,E',k)B{E',-k)) 

+t{E,k)a^{k) +r{E,k)a{-k) , (5.64) 

which obviously means that we are assuming that b^{E,k) is a linear combi- 
nation of B^{E, k), B{E, k); a^(k) and a(k). This problem was solved in [5]. 
The coefficients T{E,E',k), R{E,E',k), t{E,k) and r{E,k) depend on the 
form factor f{q) in fl5.59p and can be written in terms of the Green function 



GiE,k) = ^^^ (5.65) 

with 



CO- - E^ -U{E,k) ' 



n(^,k) = I dE'2E'0^ (5.66) 
and Eq = [Am^ + k^)^/^. Therefore, G{E, k) depends on 

through n(i?, k), since p{E, k) is given by equation ( ]5.67p . Note that p{E, k) 
depends on a{y) which is the Fourier transform of the form factor f{q). 

From f l5.67p . we see that we cannot choose arbitrarily the form factor 
/(g). For instance, if we fix f{q) = 1, its Fourier transform a{y) is a Dirac 
delta that cannot be squared thus making fl5.67p meaningless. To be in the 
safe side, we may choose the form factor /(g) to be a smooth (i.e., Schwartz) 
function. 
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For k fixed, the function G{E, k) is a function of the complex variable E 
with a cut on the real semiaxis [Eq, oo). The boundary values of the complex 
variable function G{E, k) from above to below and from below to above are 
respectively given by G+{E,'k) and G-{E,'k). 

The coefficients T{E, E', k), R{E, E', k), t{E, k) and r{E, k) in flCTjl de- 
pend on G{E, k) and hence on the operators k) and its adjoint b{E, k). 
If we use G+(-E,k) instead of G^E,^), we are using the so called incoming 
boundary conditions. Henceforth, we shall denote by bl^{E, k) and bin{E, k) to 
the solutions of (15.641) and its adjoint equation, where we have used G+{E, k) 
instead of G(£^, k) (If instead of G+(£',k), we use G_(£',k), we obtain the 
new solution bl^^{E, k) and bout{E, k) corresponding to the outgoing boundary 
conditions). 

Now, there are two possible situations. If M < 2m, the functions G+{E, k) 
and G-{E, k) are analytic with no singularities on the upper and lower half 
plane respectively [9]. On the other hand, if M > 2m, G+(i?,k) has a pole 
at = k^ + /i^ and G^{E, k) has a pole at = k^ + /i*^, where the star 
denotes complex conjugation, and 

= fx^ — ijjT . (5.68) 

The real and positive numbers fi and F depend on the form factor f{q). 

The existence of a pole of this kind implies the presence of a metastable 
state (resonance) of the system [H]. Carrying the analogy with the previous 
model further, the metastable state appears when the particle of mass M 
can decay into the two particle system described by the bilocal field. The 
condition M > 2m is of course necessary for this process to take place. 

The equation for the complex pole fic is [9] 

uj\\.)-E^- jdE-^^^=^. (5.69) 
For small values of the coupling constant A, we find [9] 

VM2 - 4m2 ^ ^ 

In order to obtain the Gamow vectors for this resonance, we have to 
obtain first the vacuum Q being annihilated by the operators 6(ii^,k). This 
can be obtained by the vacuum |0), annihilated by i?(i?,k) and a(k), i.e., 
the initial vacuum state defined by 
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B{E, k)\0) = 0, a(k)|0) = 



(5.71) 



This vacuum is not annihilated by b{E, k), i.e., b{E, k)|0) 7^ 0. As a function 
of the energy E, 6(£', k) is analytic with a branch cut in [Eq,oo) [9\. The 
new vacuum Q is required to have the following property 

kn{E,k)n = 0, (5.72) 

where 6in(-E',k) represents the function of the boundary values of 6(£', k) 
on the cut from above to below. We obtain this new vacuum making use 
of the standard theory of Bogolyubov transformations [22]. It should be a 
superposition of states with an arbitrary number of particles of B and a types 
and is obtained from the old vacuum |0) by a Bogolubov transformation 

n = e^ |0) , (5.73) 

where V depends on the creation operators B\E, k) and a^(k) and on three 
functions depending on the energy and the three momentum. Its explicit form 
is rather complicate and not strictly necessary to understand the sequel. See 
details in P]. The procedure requires a rather cumbersome calculation that 
we prefer to omit here, see [9]. Now, let us define 

$i„(E,k) = fe^j£;,k)|n), (5.74) 

which, as a function of the complex variable E, has a pole in the analytic 
continuation from above to below. This pole is located at the point = 
(k^ + /i^ — ifiVy^^ |9]. On a neighborhood of zr, the vector ^i^iE, k) has the 
form: 

$in(-E, k) = — ^ Vhi(X) + regular part . (5.75) 

E- zr 

By construction, ^p^(k.) is the Gamow vector associated to the resonance 
pole zr. Mathematically, ^p^(k.) can be rigorously defined as a functional 
on a space of test vectors, dense in the Fock space [9J, so that the following 
properties hold: 

• The Gamow vector </'|^(k) is a eigenvector of the Hamiltonian with 
eigenvalue zr, i.e.. 
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PoV^rn(k) = ^i?,<(k). (5.76) 

It decays exponentially (in a weak sense, note that f^O^) is not in the 
Fock space), so that if t > 0, 

e-'*^>Pn(k) = e-^*^«V^g(k). (5.77) 



It is important to remark that equation (15.771) implies that the Gamow 
vector v^i^(k) decays exponentially. In fact, we can decompose zji = (k^ + 
/i^ — i/iF)^/^ in its real and imaginary part: Zr = A{k.) — iB(k). The form of 
the imaginary part is given by 

For small values of F, we obviously have using the Taylor theorem 

From fl5.79p and f l5.77p the exponential decay of the Gamow vector v^i^(k) 
follows. 

5.2.1 Construction of the rigged Hilbert space. 

Since the Gamow vector v^|^(k) decays exponentially, it cannot belong to a 
Hilbert space in which the time evolution is unitary. There is one possibility 
for ip^(k) to decay under time evolution which is to extend this time evolution 
into the dual space of a rigged Hilbert space. This dual space will include 
!f^(k). The procedure is identical to what we have done in the case of the 
standard Friedrichs model although the construction should be obviously 
more complicated. 

We begin with the following space 

^ = {n'ins)® s{R^) , (5.8O) 
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for which "H^, S and iS(M^) have been defined in the previous section right 
after (CTI]) . A typical function in (nl fl 5) (g) S{R^) is of the form ip{E, k). 
For a fixed value of the momentum k, this is a function of the energy which 
belongs to Ti^ fl S. For a fixed value of the energy E, (p{E, k) e iS(R^). 

Then, if (p{E,k) is in * = (nl n S) S{R^), we define the following 
vectors, written in bra form: 

/■°° rl^krlF, 

where \ denotes the new vacuum introduced in (15.731) . Note that the above 
integral goes from Eq to oo, which means that the energy is always positive. 
Also, its is important to point out that for fixed k, the function if{E, k) is a 
Hardy function on the lower half plane that is completely determined by its 
values on the interval [Eq, oo) [9]. 

By construction, the set of the abstract vectors of the form fl5.8ip is a 
vector space The mapping 

/: ip{E,k)^\if) (5.82) 

defines a one to one correspondence from ^ = ("H^ fliS) ®5(M^) onto We 
do not have defined a topology in $ yet, but we have an isomorphism g from 
^ onto $ and a topology on which is a tensor product of two Schwartz 
topologies [92]. Then, we use g to transport the topology from ^ to This 
is a usual procedure in the theory of locally convex spaces [92] . 

Then, we want to outline in here the construction for the rigged Hilbert 
space in which are well defined the Gamow vector ^p^(k) and its time evo- 
lution. We shall skip some technical details that go beyond the scope of the 
present work. 

First of all, we construct tensor products of the form $ $ ^ ® ^, 
etc. The algebraic tensor products give topological vector spaces that have to 
be completed with respect to some topology (which in particular assures the 
convergence of Cauchy sequences) [60]. Then, let us construct the following 
sequence: 
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*0 
*2 



= c 
= c©$ 



(5.83) 

Note that, for all values of n, C $n+i- It is also true that the iden- 
tity mapping from into $n+i is continuous. Then, let us construct the 
following space: 

oo 

Y_ := IJ . (5.84) 

fc=0 

Then, we endow y_ with the finest topology that make all the following 
identity maps continuous: 

In : ^ . (5.85) 

Roughly speaking, finest means that any other topology in that makes 
all the in continuous must have less open sets. The topology in YL is usually 
called strict inductive limit topology (GOj [921 1100] . The space YL is a nuclear 
space [60] . 

We denote by the dual space of Y_. It has the property that for all 
natural n, C © © ■ ■ ■ © ($^ © ■ ■ • © ^>^) C Y^ , where the last tensor 
product has n factors [9]. 

To complete a rigged Hilbert space, we need a Hilbert space. Let us 
start with n't © L'^{M?). This space is the completion of (H^ n 5) © S{W) 
with respect to the Hilbert space topology. It is naturally isometric to the 
completion H of $ with respect also to the Hilbert space topology. The space 
"H can be looked as the result of using functions in "H^ © L^(M'^) instead of 
functions in {Wi n 5) © ^(M^) in (ICTj) . Then, 
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'Hn = C®l-i®---®{'H®---®l-i) 

and then construct the Fock space 



C 



(5.86) 



oo 



(5.87) 



fc=0 



Then, the triplet 



r_ c ^ c r 



X 



(5.88) 



is a rigged Hilbert space. The Gamow vector v^|^(k) is a continuous functional 
on y_ and therefore, it belongs to . When applied to a vector in C or 
$ ■ • ■ it gives zero, provided that the tensor product involves two or 
more copies of On $ itself, </3|^(k) is nonzero. 

The equation e**^°y_ C Y_ is valid if and only if t > [9j. Therefore, 
by duality (see Appendix), we obtain that the time evolution operator e~^*'^° 
acts on with e'^^^^Y^ C Y^ if and only if t > 0. Thus, (ISTTT]) is correct 
if and only if t > 0. 

Also the operator Pq has the property that PoY_ C Y_ so that it can be 
extended into the dual Y^ . It has the property that [Hj. 



Here we conclude this discussion on the construction of the rigged Hilbert 
space (15.881) often called the rigged Fock space. 

5.3 Priedrichs model with virtual transitions. 

The Friedrichs model with virtual transitions is probably the simplest Friedrichs 
model for unstable interaction that uses second quantization techniques. This 
is an approach which is completely different from the models discussed so far 
as we shall see in the sequel. 

5.3.1 The second quantization of the Friedrichs model. 

As we have see in section 2, the total Hamiltonian in the Friedrichs model 
has three terms H = Ho + He + Hf. 



(5.89) 
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i. ) A continuous part, He, represented by the integral in f l2.9p . This in- 
tegral is like a "linear combination" of the eigenkets, of the free Hamil- 
tonian Hq whose eigenvalues are in the continuous spectrum of Hq, which 
is the positive semiaxis = [0, cxo). Thus, the eigenkets \uj) fulfill the 
relation Hq\uj) = uj\u). For \uj), we can define creation, b}^, and annihila- 
tion operators, b^, so that from a hypothetical vacuum |0), blj\0) = \uj) and 
buj\uj) = |0). Clearly, bl^b^ are the number operators for states \uj). Obvi- 
ously, bljbuj\0){0\ = \uj){ijj\, so that we can write the integral in (12 .Op in terms 
of these operators as 

/•oo /"OO 

Hc= duu\uj){uj\ = / duublbjO){0\. (5.90) 
Jo Jo 

Then, we can omit the dyad |0)(0| in (I5.90p and simply write 

POD 

Hc= duublb^. (5.91) 
Jo 

ii. ) A discrete part, Hd, represented by the first term in (12. 9p . This 
discrete part consists in a single eigenvalue imbedded in the continuous spec- 
trum of Ho. This eigenvalue has an eigenvector that we call |1). If and a 
are the respective creation and annihilation operators for the "discrete" 
part 

HD = uJo\'i-){M=^oa^a\0){0\, (5.92) 

or simply, 

Hn = uq a^a . (5.93) 

Note that uq > 0. 

iii. ) The interaction between the discrete and the continuous parts. Hi, 
depends on both the coupling constant A and the form factor /(w). the 
explicit form is given in (I2.10p (except for the multiplication by A). In terms 
of the creation and annihilation operators following the notation as used in 
(EUD and (K93\\ . Hi is given by 

POO 

Hi = X dujf{u){a)b^ + abl). (5.94) 
Jo 

These creation and annihilation operators satisfy commutation relations 
corresponding to bosons, so that 
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[a,a^] = / ; [b^,bl] = 6{uj - u') , (5.95) 

so that the Friedrichs model can be seen as a boson interacting with a boson 
field. 

In the sequel, we want to propose a slight modification of the above model 
with total hamiltonian given by H := Hq + AV^ with 

POO 

Hq := £0 + 1^0 o-^o- + / duubljb^ , (5.96) 
Jo 

where eg > is the energy of the vacuum. The potential V is given by 

POO 

V:= duf{co)ibl + b^){a^ + a). (5.97) 
Jo 

We have chosen a real form factor which should be, in addition, square inte- 
grable: 

POO 

/ duf{u)< 00. (5.98) 
Jo 

Note that H has a continuous spectrum given by [eq, 00) and a discrete 
spectrum, which has an isolated point at Eq + ujq, which is embedded in the 
continuum. 

In order to obtain a compact form for H = Hq + XV, we need to di- 
agonalize this operator. When we use this second quantized version of the 
Friedrichs model it is useful to do it by solving the eigenvalue problem of the 
Liouville-von Neumann operator L = [H, ■] acting in the algebra of observ- 
ables generated by the operators Oi, a\, b^^ and 6^. This eigenvalue problem 
can be posed as 

LBl = ujBl ■ LB^ = -ujB^, (5.99) 

where and B^^ are functions of the previously given creation and annihi- 
lation operators respectively [HIIE]. The final result will be the diagonalized 
total Hamiltonian which has the following form 

POO 

H = Eo+ / duiuBlB^, (5.100) 
Jo 

where Eq is a renormalized vacuum energy that will be obtained later. 
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This problem can be solved as follows: We write the operators 5^ and 
B^^ as linear combinations of a^, a, h]^ and 6^ with unknown coefficients that 
are unknown functions and constants. The procedure has been explained in 
[9]. These linear combinations are then inserted in (15 .99 p . It results a set of 
integral equations for the unknown functions and constants. The solution of 
these integral equations give the mentioned coefficients. The final solution is 



bl + 2LOo\f{co) G^iu) 



oo 



bl 



h,.,i 



X < / duj'X f{uj') 

/o \uj' — (u ± w' + a; 



{oj — (jJq) + {(jJ — ojq) a 



(5.101) 



X < ^ duo'X f{u') 



(B^) in =b^ + 2uoXf{u)G^iuj) 
6,v \ {u — Uq) a"!" + (w — uq) a 



(to =F iO) to' + UJ 



2^0 



(5.102) 



Two sets of solutions called incoming and outgoing are determined by the 
boundary values of the Green function G{z) on the real axis. The boundary 
values from above to below are denoted as G~^{u) = G{uj + iO) and those 
form below to above as G~{u) = G{u — iO). The Green function G{z) for 
complex values of z is given by 



G{z) 



(5.103) 



If, in addition, the form factor f{uj) satisfies the following condition 

(5.104) 



- / du ^ > 

^ 



then G{z) is analytic in the whole plane without any singularities except for 
a branch cut coinciding with the real semiaxis [eo, C)o). Nevertheless, G{z) 
admits analytic continuations through the cut as mentioned before. 
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The new operators satisfy simple commutation relations: 



{BlU = [(5^)out, (5i,)out] = S{u - u') . (5.105) 

The next step is to write the original field operators in terms of the new 
ones. As we are mainly interested in the decay process we are in the time 
sector t > 0. As the "in" operators are ultimately responsible for creation 
and annihilation of the decaying sector, we are from now on assuming that we 
are using only incoming operators. Consequently, we shall drop the subscript 
"in" unless necessary. With the help of the properties of the function G{z) 
[6T]. we obtain that 



bi = Bl-2uo\f{u) / du'Xfico') 

Jo 



-D,v : D:„l 



u' — U — iO 



UJ' + Ul 



(5.106) 



b^ = B^-2uoXfiuj) / du'Xfiu' 

'o 



oj' — u + iO 



B„,r 



Bi, 



(5.107) 



at = -/ duj\f{uj)[{uj + ujo)G-{uj)Bl-{uj~uJo)G+{uj)B^], (5.108) 
Jo 



doo X f{uj) [{uj + Wo) G^{uj) B^-{uj- ioo) G-{uj) Bl] . (5.109) 



We note that there is a substantial difference between the representation 
in terms of a and 6^ and the representation in terms of (and their hermitic 
conjugates). The reason is clear: while the free Hamiltonian Hq has a bound 
state, and we need the operators a and to annihilate and create this state 
respectively, this situation does not arises for the total Hamiltonian H, which 
solely has continuous spectrum. Thus, the Bogolubov transformation that 
transfer equations in terms of operators of the type a, 6^, etc (often called bare 
operators) into Blj and B^^ [dressed operators) are of the so called improper 
type |2l]. In any case, it is possible to construct a new vacuum state \VL) 
from the vacuum state |S) annihilated by a and 6^. This new vacuum state 
will satisfy for all uj e [sq, oo). 
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B^\n) = 0. (5.110) 

This new vacuum state along the commutation relations flS.lOSp for the op- 
erators Bl and B^^ permit us to construct the Fock space suitable for the 
space of states of the new situation. Following the general approach for the 
Bogolubov transformation [2T] , the new vacuum state is related with the old 
one through the following expression 

|n):=e^|S), (5.111) 

where V is the so called dressing operator and is given by a quadratic form 
on the bare operators. The unknown constants and functions that serve as 
coefficients of this quadratic form are obtained by inserting it into (15.1111) 
and then using (I5.110p . The technique was already used in the previous 
model, see [9]. The final result is 



V= du' du" uo y f{u') f{u") r]{u') r]{u") 



X 



2{uo + Seo) u' + u" 



where feo is the vacuum energy shift. This result allows us to determine the 
vacuum energy Eq of the total Hamiltonian H as 

j5,=e„+ r<i.^£^!m:M=e„+fe„. (5.113) 

Jo ^0 + deo 

The function ri{u), or for complex values of the argument r]{z), is the 
solution of the factorization of the Green function G{z), a problem solved in 

ri{z)ri{-z) = G{z) . (5.114) 

The function r]{—z) is analytic on the complex plane with a branch cut in 
(— oo, — £o]- It obeys the following expressions: 
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'o 



u — z 



(5.115) 



and 



Jo Vi^) - z) 



(5.116) 



5.3.2 Time evolution for the bare operators. 



From equations f l5.106ll57l09p . it is possible to obtain the time evolution with 
respect to the total Hamiltonian of the bare operators. The derivation of 
next formulas appear in pT] and the final result is given by (&J;(0) = 6jj) 



Acot 



/•oo 

e-bl{0)-2cuoXfiu) / dcu'Xfico') 

Jo 



hi 



u' — u — iO 



[g(u,t)-gi-uj',t)]- 



b,.,' 



d 



d 



g{u,t) (5.117) 



and 



a\t) 



~ / ^^^^ 



+ / dujX^f{uj)\G{uj)\^{[{uj + uofe'^'-{uj-uofe-'^']a^ 

'0 



+ (a;'-a;o)'[e*" -e-n«}, (5-118) 



where 



2 1^ , ^ r^/^\\ Jzt 



g{z,t):= I duj2uj,X^ f{uj)\G{uj)\ ^ 

Iq \uj — z uj + z 



+ G{z) e 



(5.119) 
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One can prove that g{z,t) is an analytic function of z except for the 
infinity, see [HI]. In terms of t, it has the following asymptotic behavior: 

g{uj, t) I — > -e-^^* C^iu) as t ^ ±00 . (5.120) 

Also, we can obtain [HI] the asymptotic behavior, in terms of t, of the bare 
operators: 

bl I — > (5i)out as t^±oo, (5.121) 

in 

and 

a\t)i — ^0 as t^±oo. (5.122) 

Equations 05.1211) and fl5.122p are valid in a weak sense so that no con- 
clusion about the products a'^(t)a{t) and a(t)a^(t) can be made. Therefore, 
there is no contradiction between fl5.122p and the commutation relation 

[a{t),a^{t)] = l, (5.123) 

valid for any finite value of t. 



5.3.3 Resonance poles. 

In order to find resonance poles, we have to go back to fl5.103p . In fact, 
G{z) (15.1031) has the same structure as ri~^{z) in (12.141) corresponding to the 
standard Friedrichs model. Therefore, the analytic continuation of G{z) has 
two complex conjugate resonance poles located at the points zr and z^. In 
the first order of approximation in A, we find that 

ZR^Uo + i'^, (5.124) 

where 

7 = 27rAV'(c^o) (5.125) 

and 

poo \2 f2(^\ 

Uo-PV I duj !^ , (5.126) 



- ^5 



where PV stands for principal value. 
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5.3.4 Why Priedrichs model with virtual transitions? 



As an example of the application of the solution of this second quantized 
Friedrichs model, we shall consider the time evolution of the number of bare 
photons. Starting with the bare vacuum |S) as the initial condition, the goal 
is to calculate the quantity given by 

n^t) := {E\blb^\E) , (5.127) 

where obviously blj b^^ is the number operator for photons with energy u. 
Then, using (15.1 17p and its Hermitian conjugate, we obtain the following 
result: 



Jo + w ) 



«^ - Wo ) g{-uj,t) 



. (5.128) 



This expression is exact and it is valid for any value of A verifying fl5.104p . 
In the weak coupling case (small A and constant tA^), we can obtain the 
following formula [611 182] : 



n^(t) = 2 7^^^^^ [1 - cos((u; + uo)t)] + O(A^) 



(5.129) 



{UJ + LOq^ 

This formula reflects the initial grow of the number of bare photons due 
to virtual process. According a point of view [82], this corresponds to the 
formation of a photon cloud around an atom. Also, the first term in (15.1291) 
gives the probability of emission of a virtual photon. 



6 Other types of Friedrichs models. 

New kinds of Friedrichs models have been studied in the literature. Along 
the present section, we discuss two more versions. In the former, we have 
a discrete Friedrichs model, which has been used to approximate systems 
with a continuous spectrum by systems with a discrete spectrum [53] . In the 
other, the internal channel has already continuous spectrum and this allows 
to the appearance of a resonant branch cut in the analytic continuation of 
the reduced resolvent [HI], as we shall see. 
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6.1 A discrete Friedrichs model. 



We begin with a harmonic oscillator, the system (in the language of thermo- 
dynamics) or the internal channel (in the language of the Friedrichs model), 
interacting with a bath of harmonic oscillators (the external channel), for 
which the free Hamiltonian is given by 

n=l ^ ' 

Following [53], we use capital letters for the variables in the system and 
small letters for the elements of the bath. This makes a proper notational 
distinction. We need to define an interaction between both system and bath 
(or internal and external channels), and the proposed interaction reads [53] : 

i/,:=|:c„(.Y.,. + ^^). (6.2) 

where the c„ are real (due to the need for hermiticity of Hi these variables 
must be real) coupling constants. Note that Hj is non local. The total 
Hamiltonian is H := Hq+Hj. If we define the usual creation and annihilation 
operators as (we omit the constant h or make it equal to one): 



2 V 2mnUJn 



^n + l\h, Pn, (6.4) 



which satisfies the canonical commutation relations given by 

[B,B^ = I ; [bn,bl]=Snm, (6.5) 

where Snm is the Kronecker delta and this equation is valid for any 71,771 = 
1,2 . . . , N. All other commutators vanish. In terms of these operators, the 
total Hamiltonian H has the following form: 

H = n {^B^B + -j+J2un [bib + 2 J + E 9niBbi + B^bn) , (6.6) 
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where gn '■= Cn/\/ MQmnOOn- This kind of couphng has been previously 
considered, for instance in [69] . 

Let us assume that the vacuum state is defined as 

|0) := |0)® |0,...,0), 

where the first term |0) in the above tensor product is the ground state for the 
system and |0, . . . , 0) is the ground state of the bath, in which we distinguish 
each harmonic oscillator. For instance |ai, 02, ... , a^) is the state of the bath 
with the first oscillator in the state Oi, the second in the state 02, etc, up 
to the latter in the state a^v- In the present case, ai = 02 = • • • = ajv = 0. 
Define now the following vectors: 

\n) = El'lO) = |1) |0,...,0) (6.7) 

\^k) = 6l|0) = |0)®|0,...,l,...,0), (6.8) 

where the 1 in (16.81) is placed at the k—th site. With this definition, the total 
Hamiltonian becomes: 

N N 

H = +J2^n\0On){0On\ + ^ f7„(|fi)K| + \Un) {^1) +C, (6.9) 

n=l n=l 

where C = Q/2 + E„=i(wn/2) is a constant!^. The Hamiltonian in (16. 9 p has 
already the structure of the Hamiltonian of the Friedrichs model (except for 
the constant that we may eventually drop). In fact, (16. 9p has a free Hamilto- 
nian given by + '^n=i^n\(^n){(^n\ plus an interaction that connects 
the internal and external channels given by J2n=i fi'n(|^)(i^n.| + |i^n)(^|)- The 
contants Qn play the role of form factors. In addition, we impose the obvious 
condition that fl should not be identical to none of the Uk- 

In order to solve an eigenvalue equation, similar to (I5.99p . we first make 
the following Bogolubov transformation: 

N 

Ck := AkB + ^aknbn, (6.10) 

n=l 

^^Note that, save for the constant C, the Hamiltonian in (|6.9I) is given by i/|0)(0|, where 
H is the Hamiltonian H given in ()6.5p . The situation is quite similar than those discussed 
in Section 5.3. 
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where the coefficients A/^ and akn, n = 1,2, . . . , N, k = 0,1,2, N are 
chosen in such a way that the transformation be unitary and the following 
commutation relations hold 

[cfc,c|] = 6ki , (6.11) 

and all other vanish. In addition, we require transformation fl6.10p to diago- 
nalice the total Hamiltonian H as: 

N 

H = Y,(^kclck + C. (6.12) 

k=0 

This requirements completely determine the transformation f l6.10p . In fact, 
if \ak) := cj.|0), the coefficients in f l6.1Up are given by 

Ak = {ak\^) ; ttkn = {ak\(^n) ■ (6.13) 

The vectors la^) are orthonormal. We want to obtain the eigenvalues of the 
renormalized Hamiltonian := H — C, C being the constant that appears 
in the expression of H in fl6.9|6.12p . 

The operators Ck, B and do not commute with the total Hamiltonian 
H (or and their commutators give: 

TV 

[ck,H] = akCk = akAkB + ^akaknbn, (6.14) 

n=l 

TV 

[B,H] = QB + Y,9nbn, (6.15) 

n=l 

[bn,H] = uJjn + gnB. (6.16) 

These equations give: 

TV 

(6.17) 

n=l 

akttkn = AkQn + UJnakn- (6.18) 
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From f l6.18p one obtains 



an - 

If we carry (16.191) into 06.181) . we obtain 



akn = • (6.19) 



N 2 

Akiak-n-Y,^^] =0, (6.20) 



n=l - ^« 



which yields 



N 2 

ak-n-y" — — — = 0, k = 0,l,...,N. (6.21) 

Equation f l6.2ip is transcendental and can be solved by numerical meth- 
ods. 

As in any Friedrichs type model, we are interested in the reduced resol- 
vent. Let P be the orthogonal projection into the external channel, the one 
dimensional vector space spanned by the vector and let Q := I — P, 
which is the projection into the N dimensional vector space spanned by the 
{l^n)}, n = 1, 2, . . . , A^. In particular, P = The reduced resolvent is 

given by the operator 

P-l—P = ^-\z)P. (6.22) 

As in the ordinary Friedrichs model, P projects into a one dimensional sub- 
space so that the identity in f l6.22p holds. In order to find the function ri~^{z), 
we shall make use of a formula that holds in this case and that we show in 
Appendix C (see also 138]). This is (IB.ISP which reads (the original form 
includes the coupling constant A that here is taken equal to one): 

Pj^^P=[G{z)]-'P (6.23) 



with 



G{z) = {z~ Ho)P + PV P - PV Q — \-QVP. (6.24) 

z — Hq 
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Then, we get G{z)\^l) and this immediately gives f l6.23p as we shall see. Note 
that the potential V is here given by = J2n=i 9ni\^) {'^n\ + |i^n)(^|) or in 
matrix form: 



V 



( g\ ■■■ gN \ 

9i 







\ 9n 



(6.25) 



This matrix is written in the basis {\^), \uJi), ■ ■ ■ , \ujn)}- Note that 



/ 

9i 



\ 9n 



9i ■■■ 9n 







/ 1 \ 



/ \ 



91 



\9N J 



Let us apply the last term in ( I6.24p to \ fl). This gives: 



(6.26) 



PVQ^—QVP 
z- Ho 



PVQ^—Q 

z- Ho 



We know that 





/ \ 

91 

\9N J 



( \ 

91 



PVQ —QV 



z-Ho 



f 1 \ 



/ \ 



PVQ 



z-Ho 



91 



\9N J 



N 



^9n\^n) 



n=l 



\9N / 

and that Ho\ijJn) = oOn \oJn)- Then, fl6.27p equals to 



(6.27) 
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PVQ 



( \ 

Z—L>J\ 

\ I 



PV 



/ \ 

.91 

\ I 



( _gi_ \ 



p 



\ 



J 



The other terms are: 



(6.28) 



PVP 



/I \ 









[z - Ho)P 



( 1 \ 



voy 



( z-9.\ 




Therefore, 



N 2 



n=l 



2; — CJn 



and then, 



[G{z)]-^\^)=[z-^-Y, 



N 9 



-1 



n=l 



2; - uj„. 



Thus, obviously, 



Af 2 

77(2) = ^ - - V 



(6.29) 



(6.30) 



(6.31) 



(6.32) 



Compare fl6.32p with fl2.14p . We know that the equation rj{z) = has real 
solutions given by ao, . . . , a^v. 

This model has been used to compare the results given by a system with 
discrete spectrum and another similar with continuous spectrum. The usual 
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argument consists in taking the limit of the sum in f l6.32p in a way that the 
sum is transformed into an integral of the type that it appears in fl2.14p . 
Some results are discussed in |53l [981 ES] . 

6.2 A generalized Friedrichs model with branch cuts. 

In this section, we want to present a generalized form of the Friedrichs model 
in which, instead a single or multiple resonance pole, a continuous set of 
singular points, a branch cut, appears. This construction may appear rather 
artificial, although we intend to write an explicit example of a continuous 
resonance and a method to produce it. This example has already appeared 
some time ago pi] . 

In the simplest form of the Friedrichs model, we have a free Hamiltonian 
with one eigenstate plus a continuum. Therefore, the model has two channels, 
one infinite dimensional corresponding to the continuous spectrum of Hq, 
the external channel, and the other with dimension one corresponding to the 
stationary state of Hq, the internal channel. The mutual interaction of these 
two channels produces the resonance. This pattern is shown in any kind of 
generalization of the Friedrichs model. 

In our case, both channels are infinite dimensional. The external chan- 
nel has also a continuous spectrum coinciding with the positive semiaxis, 
although now this continuous spectrum has an homogeneous infinite multi- 
plicity (or degeneracy). The internal channel is here a I2 type infinite di- 
mensional Hilbert space. Both are believed to be minimal requirements to 
produce a continuous resonance. 

Thus, the Hilbert space is divided into two mutually orthogonal parts: 



Needless to say that the terms "in" and "ex" stand for internal and external 
channels respectively. The Hilbert space for the internal channel is defined 

as 



Here, hC^) is a space of sequences of complex numbers labelled by the set 
of the integers with the following condition: if {a„} is one such sequence, it 
must satisfy: 



ex 



(6.33) 




(6.34) 
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|a„P < oo (6.35) 

n.=— oo 

li A = {an} and B = {6„} are two sequences with the property fl6.35p . we 
can define the following inner product of them: 

oo 

{B\A)= J2 Kan. (6.36) 



With the scalar product fl6.36p . hC^) is a Hilbert space 

In order to construct the Hilbert space corresponding to the external 
channel, let be, for each integer value of n, a copy of the Hilbert space 
L2(M), so that Xn = L^(M) for all = 0, ±1, ±2, . . . . Then, let us define 

oo 

-^ex .= x„ . (6.37) 

n— oo 

After this definition, the elements of T-L'^^ are sequences {fn{x)} of square 
integrable functions on the real line with the property: 

oo 

E ll/n(a:)|P<oo. (6.38) 

n=— oo 

The norm in fl6.38p is the well known L^(]R) norm. The scalar product in 
'H'^'^ of two such sequences, f = {/n(x)} and g = {gn{x)} is given by 

oo «QQ 

(f|g):= J2 / f:{x)9n{x)dx. (6.39) 

n=-oo 

We see that the external channel is multiple valued. Alternatively, it can 
be looked as a superposition of infinite simply valued external channels [91] . 
It is just a matter of language. 

As for the Hilbert space, the free Hamiltonian can be split into two contri- 
butions on the internal, if™, and on the external channel, H'^^, respectively, 
so that 

Ho = H'"' ® . (6.40) 

Let us define first. If 
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^:={iJn}el2{Z)=n'\ (6.41) 
Hip will be the sequence in /2(^) such that its n— th term is 

:= a,Pn - ipn+i - ^n-i , (6.42) 

where a is a fixed real number. We can prove that if™ is a bounded (con- 
tinuous) operator on /2(Z) and that is has a doubly degenerate continuous 
spectrum given by the interval [a — 2, a + 2] [91j . 

On the other hand, H^'^ must be an unbounded operator on "H*^^, as we 
requiere that its spectrum is = [0, oo). For unbounded operators, we 
need to define the domain of the operator (the space in which the operator 
acts) plus its action on the domain. For W^^, the domain is defined as 

oo 

= , (6.43) 

n=— oo 

where, for all n G {. . . , /c, . . . , —2, —1, 0, 1, 2, . . . , fc, . . . } = Z, we have 

ir„ = iy|(M), VnGZ. (6.44) 

The Sobolev space Vr2^(M) is the space of twice differentiable functions 
f{x) from the real line R into the set of complex numbers C such that 

'{|/(x)p + |/"(x)r}rfx<oo. (6.45) 

This Sobolev space 14^2^ (M) with the square norm given by (16.451) is a Hilbert 
space. In addition PF2^(]R)is dense in L^(]R) (being given a function g{x) G 
L^(]R) and an e > 0, there is always an g{x) G W2 such that ||/ — < e 
in the norm), which implies that the domain of H^^, ^{H^^) is dense in 

Once the domain of has been defined, we define the operator H^^ by 
giving its action on an arbitrary f G T>{H^^). If f = {/n(x)}, we have 
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It is a technical exercise to show that H*^^ is essentially self adjoinll^ 
and that it has an absolutely continuous spectrum given by [0, oo) with an 
homogeneous infinite degeneracy!^. 

As a consequence, the Hamiltonian Ho = if™ + H^^ is self adjoint with 
domain © V{H^^) and a spectrum given by 

a{Ho) = U = [a - 2, a + 2] U [0, oo) , (6.47) 

with infinite degeneracy. 

This Hq is the unperturbed Hamiltonian. We need to find a perturbation 
such that the part of a (if™) that intersects a{H'^^) = [0, oo), which obviously 
is [0,a + 2], is "disolved" and becomes a resonance in form of a branch cut 
in a suitable analytical continuation of the resolvent. 

To find this perturbation, we shall use a technique developed in [SSI 1^ . 
It has three steps: i.) Finding a dense domain P C V{Ho) such that the 
restriction, Hq, of iio into V admits infinite self-adjoint extension^; ii.) We 
find the adjoint Hq of Hq. This adjoint does exist as Hq is densely defined 
(its domain is dense in "H) and this adjoint is an extension of Hq (for any 
/ in the domain of Hq , HqJ = Hq f; iii.) Providing suitably boundary 
conditions on the domain of Hq\ we find an extension of Hq of the form 
Hq + 'yV (which is also a restriction of Hq , i.e., it acts as Hq on a domain 
smaller than the domain of Hq^). This operator H^ := Hq + ^V is self adjoint 
and V is the desired potential. Note that H^ has the typical form of the total 
Hamiltonian in a Friedrichs like model. 

In order to fix Hq it is enough fixing its domain V{Hq), since Hq is a 
restriction of ifo (and therefore iiof = HQi, for all f G T>{H'q)). This domain 
is giving by 

oo 

V{H'q):=W®D- i):=0i)„, (6.48) 

n=— oo 

where Dn ■= C^(M), for all n G Z, so that all the spaces Dn are copies of the 
same C^(M). This C^(M) is the vector space of all infinitely differentiable 
functions that vanish on a certain neighborhood of the origin, which does 

^■'This means that iJ'" has one and only one self- adjoint extension and that therefore is 
a well defined quantum observable. 

''These are technical concepts that we do not want to discuss in the paper. See [20 11101) . 

^^In the mathematical language, Hq must have equal deficiency indices different from 
zero. See Appendix B. 
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not need to be the same for all functions in C^(M). One can prove that 
V{Hq) is dense in the total Hilbert space "H and therefore, the adjoint H'q of 
H'q is well defined, although H'q is not essentially self adjoint because it has 
both deficiency indices equal to cxd. In order to prove this latter statement, 
it seems convenient to finding Hq first. Its domain is given by 

oo 

V{H'^^)=V}''®C- C:=0C„ (6.49) 

n.=— oo 

with 

Cn ■■= W^|(M/{0}) ; G Z . (6.50) 

The functions in the space iy2^(M/{0}) are like the functions in W21 now 
allowing our functions and their first derivatives to have a finite jump in the 
origin. The norm square in VF2^(]R/{0}) is exactly as in fl6.45p and obviously 
W2 C Vr2^(M/{0}). Thus, f is in C if and only if f is a sequence of the form 
{fn{,x)} with fn{x) G C for all n G Z. 

Once we know the domain of Hq , we can define Hq itself. For any 
^Ij G W"", H'q^iIj = W'^ip = Hqi/j. For any f = {fn{x)} G C, 



{(i/-f)„} := ^ } , (6.51) 



exactly as in f l6.46p . Observe that Hq is an extension of Hq. In addition, 
one can show that there exists a relation between the domain of Hq and the 
domain of its adjoint, which is given by 



V{H'o^) = V{H'^)+N++U. , (6.52) 

where: i.) + means direct algebraic sum, in general nonorthogonal and ii.) 
A/± are the deficiency indices of H'q. Both are infinite dimensional and thus 
equal and therefore, H'q admits an infinite number of self adjoint extensions. 
The form of N± is obtained and discussed in [91]. 

Note that an arbitrary vector $ in H = l-L™ © can be written in a 
column vector form as: 



$ = = . (6.53) 
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Now let us assume that ^> G V{Hq). In this case, ip is arbitrary in 
"H™ = ^2(^) and f G C. As in the case of the ordinary Friedrichs model, the 
perturbation will put in interaction the external with the internal channels. 
As in the formalism given in [661 EZ], we need to impose certain type of 
boundary condition^ to obtain H^. These boundary conditions, for each 
n G Z, are the following: 



/n(0+) - /n(O-) = -7<^„ (6.54) 

/n(0+) = fLiO-). (6.55) 

In flCT]) and fl635D . /„(0±) and /^(Oi) are the limits to the left (with +) 
and to the right (with — ) of the functions /„(x) and its first derivative fn{x) 
for all n G Z. These limits must exist, but they could be different as each 
Uix) G Vr|(M/{0}). 

For each real 7, these boundary conditions define one self-adjoint exten- 
sion of Hq that we call with domain: 

V{H^) := C V{H'q^) fulfilling (EMD and ■ (6.56) 

On $ G V{H^), the action of is given by 




(6.57) 



with 



^„ = (i/>)n + 7/:(0) (6.58) 

fnix) = -C{x), (6.59) 

for each n G Z. Here, ifn and fn{x) are the n— th component of ^ and f 
respectively. 

^^As a matter of fact, in order to obtain the self adjoint extensions of a symmetric non 
self-adjoint operator with equal deficiency indices, we usually do it by imposing boundary 
conditions 87 . What is characteristic in this theory is that the extension gives a potential 
(in general singular) to be added to the original unperturbed Hamiltonian Hq 
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Each of the is certainly self-adjoint [HI]. To give meaning to as 
Hq + 7V^, let us consider $ G V{H^) and S G V{Hq) and the following 
difference: 

(iJ-,$,H)-($,/7o3). (6.60) 

Obviously, as is self adjoint, if^ can be written in the form Hq + '-jV if 
and only if (16.601) is equal to 

7(1/$,H). (6.61) 

If 

{in} 

g / V {9n{x)} 

where ^ G /2(^) and g G -D, one gets [91] 



(i/^*,H)-($,ifoH)=7 E /r(0)a + 7 E <C(0). (6.62) 

n=— oo n=—oo 

Note that the boundary condition f l6.59p gives a unique value for each /^(O) 
and that ^„(0) for each n G Z is well defined because ^n(a^) £ We 
can obtain an explicit for of V using the following matrix kernel form: the 
difference {H^^, H) — HqS) is a sum in M and n of terms of the form: 

/ Vmn Vmniy) \ f in \ 

«,/;(!/)) =/r(o)en + <^;(o). 

/ \ 9n{x) ) 

(6.63) 

This yields: 



(6.64) 



Thus, the construction of the generalized Friedrichs model is finished and 
it remains to find its resonance behavior. This will be done in the next 
subsection. 
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6.2.1 The resonance as a branch cut. 



The procedure to find the resonance branch cut foUows exactly the same 
procedure as in the case of the simplest Friedrichs model. In the latter, we 
studied the singularities of the analytic continuation of the complex func- 
tion given by (12.12112.131) . which expectation values of the operator reduced 
resolvent (12. lip . In this form of the Friedrichs model, the internal channel 
is the one dimensional Hilbert space spanned by the vector |1). Thus, the 
generalization of (12. lip to our case is clear: here P must be the projection 
into the internal channel "H™. Then, the reduced resolvent is now given by 
P{H^ - zI)-^P. 

To carry out our analysis, we shall use the kernel form of the resolvent as 
follows: 



{H, - zl)-' I I = I M = f I . (6.65) 




I oo 

1 n=— oo 



[ oo 

— oo 



{9n{x)]Z 

li^ = {<^n};r=-oo e km = W'' and f = {/„(x)}- _^ G we write: 

<^ /"OO 

in = E ^-^-+ E / R^^n{y)fm{y)dy (6.66) 



gn[X) 



oo oo „oo 

E Rlniix)'Pn+ E / Rl^iix,y)fm{y)dy. (6.67) 



m=— oo m=— oo 



The kernels Rmn, R^mLiy)-, Rmn{x) and Rmn{x,y) all depend on the com- 
plex variable z, although we must stress that the variables of these kernels, 
here denoted as x and y, do no have anything to do with the real and imag- 
inary part of z. 

For every /c G Z, we define the following sequences: 

Rk-. = {R'nhn=-oo (6.68) 

^kix): = {/?2(x)}jr=_oo (6.69) 
and the following vector: 
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Rk \ 

(6.70) 

Rfc(x) / 

It is not difficult to show that, for each G Z, S^. is in the domain of 
[9T] . For our purposes, it is sufficient to consider k = and this gives the 
vector Hq. For simplicity, we can write 



_ .— ^0 




R 
R(x) 



(6.71) 



Let 6 be the sequence {^no}5JL-oo5 where 6no is the Kronecker delta. This 
sequence has all its entries equal to zero except for 6oo = 1. Then, we can 
prove that 



(H, - zI)S =\ I , (6.72) 

which yields: 




iW''R)n + -fKiO)-zR„ = Sno, VnGZ, (6.73) 

-R'^{x) - zR^{x) = 0, VnGZ. (6.74) 

We first solve the differential equation given by f l6.74p and find a solution 
in l^2^(M/{0}) with the boundary conditions (I6.54f6.55p . This solution, with 
k = y/z, is 

r Cn 6*'=^ ; X > 
Rn{x) = < . (6.75) 

[ C; e"^^^ ; X < 

The boundary condition ( I6.55P reads i?^(0+) = -R^(0— ) in this case. It 
implies that C„ = C^. Condition (16.56^ is here 2C„ = — 7-Rn- This gives: 

K(0) = ikC^ = --^Rn^fz. (6.76) 
If we carry (16.761) into (16.731) . we have that 



80 



[{W- - ^)R]„ - '-^""Rn^ = <5„o , (6.77) 
which after the definition of if" given in f l6.42p , one gets 

|a - Z - ^7^V^| Rn - Rn-l - Rn+l = ^nO • (6.78) 

The identity f l6.78p is vahd for any integer n and has the following solution: 



jiu\n\ 

Rn — — 

where 



Rn = TT— ' (6-79) 
2t sm u 



2 cos z/ = a — ^ — -^7^1/^ . (6.80) 

Observe that the equation ( 16.80p has complex solutions in the indeterminate 
u. The dependence on the variable z is contained in u = z/(z). 

Once we have obtained the solution for /c = in f l6.70p . it is immediate 
to obtain the solution for arbitrary k and one gets 

piu\n—k\ 

Rnk = 7T— ■ (6-81) 
2t sm u 

Let us go back to the reduced resolvent written in the operator form 
P{H^ — z)~^P. Since P is the projection onto the space "H™, the operator 
pIh^ - z)-^P acts only on W''. Then, if ip := {Vn}n=~oc ^ we have the 
following expression for the n— th component of P{H^ — z)~^Pip, which is in 
H'" if z^ (r{H^y. 

oo 

[P{H^-Z)-'P^]^= J2 Rmn{z)¥>m. (6.82) 



m,=— oo 



Let us call 6k to the sequence {5mfc}m=-oo' where 6mk is the Kronecker delta. 
This sequence has all its terms equal to zero except the k—th which equals 
to one. In this case, (I6.82p should be written as: 

oo 

[P{H^-Zy'P6k]n= Yl Rmn{z)6^k = Rnkiz) , (6.83) 



m=— oo 
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where Rnk{z) is as in fl6.8ip . Formula fl6.83p can be written in compact form 
as: 

P{H.,- z)-^P5k = B.k] ykeZ. (6.84) 

The set {6k}, /c G Z of sequences in "H™ = hC^) form a orthonormal 
basis in /2(Z) [12]. Then, as P{H^ — z)~^P is a bounded operator on "H™ if 
z ^ a{H^), the reduced resolvent P{H^ — z)^^P is completely determined by 
its action on the sequences {5fc}, for all /c G Z. For A; = 0, this formula gave 
fl6.77p . Same manipulations give for any A; G Z the following equation: 

^i/in _ ^ _ ^7275^ Rfc = ; V A; G Z , (6.85) 

or equivalently: 

Rfc = (^i/'"' -z- ^7^v^^ ^fc ; V A; G Z . (6.86) 

If we combine (16.860 with (16.840 . we find that 

P{H^ - z)-'P = (^/f'" -z- ■ (6.87) 

We know that (iJ™ — zl)~^ ia analytic at all points except for those in the 
spectrum of ff™ which coincides with [a — 2, a + 2]. Then,if 

giz) ■.= z^]p'^z, (6.88) 
the operator function given by P{H^ — z)^^P is analytic except for 

g-\a-2,a + 2]. (6.89) 

Therefore, the points for which the partial resolvent P{H^ — z)~^P is not 
analytic are given by the equation 

t = z + ^J^^/^■, Vt G [a-2,a + 2]. (6.90) 

To study formula f l6.90p . it is better to come back to the momentum 
representation through the transformation k = ^Jz. Then, (I6.90p can be 
written as 
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\/t e [a -2, a + 2] 



(6.91) 



+ -r 

2 ' 



which gives 



k±{t) 




VtG [a -2, a + 2]. 



(6.92) 



If the spectrum of if™, a{W^ = [a — 2, a + 2], has a negative part given by 
[ti, ^2], where ti = a— 2 and ^2 = a+2 if a+2 < and otherwise, the possible 
values of z = fc^ are real, as k is purely imaginary, and belong to the inter- 
val (ti), fc^(t2)]- Also, those z belonging to the intervals (ti), A;^(ti)] 
and [^^(^2)? ^+(^2)] have multiplicity one and those to m[/c^ (ti), k'^(t2)] have 
multiplicity equal to two. 

The values of t in the interval [to, a + 2], where to = max{7^/16, a — 2} 
give a pair of branch cuts, as we can deduce from (16.921) . These are resonant 
branch cuts and this expression can be understood in the following way: The 
internal block of the resolvent (ii™ — zl)~^ is an analytic bounded operator 
function on the complex plane in which the spectral cut at [a — 2,a + 2] is 
removed. The total resolvent has no bounded operator limit on the cut ^91] 
and, hence, it cannot be extended through the cut in general. However, the 
internal block of the resolvent has a bounded limit on this cut as Im z > and 
can be continued from above to below. This continuation has singularities. 
In this model, one finds instead the usual resonance poles, branch cuts. On 
these cuts, the internal block of the resolvent becomes unbounded. These 
cuts are associated, by analogy to the standard Friedrichs model, to resonance 
phenomena. 
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A Rigged Hilbert Spaces 



Following the spirit of [15j, we are going to present the Gamow vectors for 
the Friedrichs model as functional in a rigged Hilbert space (RHS). By RHS 
we mean the following: 

Let T-i be an infinite dimensional Hilbert space. Let $ be a subspace of 
T-i with the following properties: 

1. - The subspace $ is dense in 1-L. This means that for any vector 

and any neighborhood V of ip there exists & ^ in V. In other words, any 
vector in H can be approximated, with arbitrary accuracy, by another vector 
in $ with the topology on T-L. 

2. - The subspace $ has its own topology, r^, (compatible with the struc- 
ture of vector space), which is finer than the topology that $ inherits as 
subspace of Ti (we shall denote by the Hilbert space topology on Ti). 
This means that with r^, $ has more open sets than with r-^. Equivalently, 
any convergent sequence by r$ is convergent under r-^, but the converse is 
not true. 

3. - Once the space $ has been determined, let us consider the set of 
mappings F from $ onto the field of complex numbers C, with the following 
properties: 

i. ) Each of the F : $ i — > C is antilinear: For any and any pair 
of complex numbers we have: 

F{ai: + f3^) = a*F{ij) + /3*F{^), (A.l) 

where the star denotes complex conjugation. 

ii. ) The mapping F is continuous. In particular, this means that if {^Pn} 
is a sequence of vectors in $ which converges with r$ to the vector ip E ^, 
then: 

F{ip^) ^ F{ip) , 

in the field of complex numbers 

iii. ) We shall call Junctionals to the set of mappings with the above 
properties. The set of functionals is a vector space, provided that we define 

^^This conditions is necessary although not sufficient for the continuity of F. It may 
happen that $ be a nonmetrizable space, i.e., a space in which the topology cannot be 
given by a distance. Also this condition suggest that, at least, $ should be sequentially 
complete in the sense that any Cauchy sequence have a limit in 
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its sum and product by complex numbers as follows: 

{aF + f3G){ij) := aF{iP) + /3G'(^) , 

for all G Here, F and G are functionals and a and /3 arbitrary complex 
numbers. The vector space of functionals is denoted as 

iv.) For convenience, the space is usually endowed with a topology 
(weak, strong, Mackey), a fact that we shall not discuss her. See and 
references quoted therein. 

V.) Let be an arbitrary vector in the Hilbert space Then, define the 
mapping F^ as follows: 

F^(^) := (^1^) , (A.2) 

where ip is an arbitrary vector in $ and ( — |— ) denotes the scalar product in 
l-L. One can show that F^p is a functional in (the antilinearity is in fact 
straightforward). Furthermore, if 7^ 0, then F^ 7^ as a consequence of 
the Riesz theorem |T9l|3T]. Then, we shall identify the functional F^ and the 
vector ip for notational convenience. 

Observe that, after this identification, the space of functionals con- 
tains the Hilbert space "H, so that we have the following relation between 
these three spaces: 

$ C C . (A.3) 

This triplet of spaces is called a rigged Hilbert space or a Gelfand triplet. The 
topology on Tx, is chosen so that, on "H, r-^ is finer than r^- The space 
$ is often called the space of test vectors or test functions (if $ is a space of 
functions). 

So far, the definition of RHS. It is interesting to comment briefly some of 
the properties of operators on RHS. 

First of all a question concerning notation. Following the usual habit of 
physicists, we shall use the Dirac notation and write for any F G ad any 

= F{^) . (A.4) 
Note that this notation is consistent with (lA.2p . 

Let if be a densely defined operator on Ti (and therefore with a well 
defined adjoint operator if^) such that 
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i. ) The space $ reduces H, i.e., for every ■?/' G Hip G $ (we can 
describe this property as C 

ii. ) The operator H is continuous on $ with r$. 

Then, if can be extended as a (weakly continuous) hnear operator into 
the dual by means of the duality formula: 



where F and if) are arbitrary in and $ respectively I. 

It is important to remark that the eigenvalue problem can be posed for 
the extension of H on the antidual In particular, Gamow vectors will 
be, in this context, eigenvectors of the extended Hamiltonian, to a certain 
antidual, with eigenvalues given by the resonance poles. 

It is also interesting to recall here a Theorem due to Gelfand and Maurin 
on the completeness of generalized eigenvectors of a self adjoint operator [STJ 
[76] . This result gives a rigorous meaning to the Dirac generalized expansion 
of a self adjoint operator |12]. It states that being given an arbitrary self 
adjoint operator on an infinite dimensional Hilbert space, a rigged Hilbert 
space can be found in which this generalized expansion exists. A more precise 
form of this theorem can be presented as follows: 

Let "H be an infinite dimensional Hilbert space and A a self adjoint oper- 
ator on "H. Without loss of generality, we can assume that the spectrum of 
A is purely continuous. Then, there exists a RHS $ C "H C such that 

i. ) The space $ reduces A. This means that for any vector G 
Aif G Equivalently, A<^ C The operator A is continuous on $ and it 
can be extended, using the duality formula (lA.Sp into the antidual The 
extension to the antidual of A, that we shall also call A, is continuous with 
respect to the (weak) topology on 

ii. ) There exists a measure /x(A) supported on the spectrum of A, cr(A), 
(it is zero outside cr(y4)), such that for almost all A G cr{A) (with respect to 
the measure /i(A)), there exists a functional |A) G such that 



We can apply this |A) to any G When A runs out the spectrum of A, 
(7(A), we obtain a function {^\\) G L'^[a{A), djJi^X)]. The mapping U : ^ i — > 

"'^^Some authors distinguish this extension with the notation [23] 



{F\Hi;) = {HF\iP), 



(A.5) 



A|A) = A|A). 



(A.6) 
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L'^[(t{A), dfi{\)] given by Uip = can be extended to be unitary from Ti 

onto L2[a(A),c//i(A)]. 

iii.) For any pair of vectors (p,ilj G the following completeness relation 
holds: 



icp,^p)= / (v.|A)(A|^)dMA), (A.7) 

where (Xlip) = (V'lA)* and the star denotes complex conjugation. 

Furthermore, if /(A) is a (measurable) bounded complex function on 
(t(A), we have in addition that for any (f,ip & ^ 



{^J{A)^)= / f{X){v\X){X\^)df,{X). (A.8) 

Ja{A) 

Formula ( lA.Sp provides the functional calculus for A. It can be also 
extended to the case that /(A) be a polynomial. In this case, f{A) is the 
corresponding polynomial in the variable A. 



A.l Examples of Rigged Hilbert spaces. 

There are several examples of RHS that are useful in Quantum Mechanics 
as well as in chaotic systems [1]. To begin with, we introduce the simplest 
and more fundamental of all, the Schwartz space. 

The Schwartz space S is the vector space of all continuously indefinitely 
differentiable (derivable to all orders) functions from the real line M into the 
set of complex numbers C with the following property: If f{x) G S, then: 



lim 

XI— f itoo 



dx^ 







for all m,n = 0,1,2, ... . 

The space S has the following properties: 

1. - Each function f{x) G 5 is square integrable. In addition, 5 is a dense 
subspace of (M) . 

2. - S has its own metrizable topology inhereted from L^(M). 

3. - The antidual space of S, contains L^(M) and, with the weak topol- 
ogy on , L^(M) is dense in 

4. - Consequently, the triplet 



'X 
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SCL' 



(A.9) 



is a RHS. 

The notion of Schwartz space can be extended to more dimensions in the 
configuration space. Thus, let us consider the vector space of the indefinitely 



differentiable functions f{xi, 



lim 

Ixl Ih- ^oo 



from M"^ into C, such that 



fix. 



0. 



(A.IO) 



dx"^ . . . dxq '• 

where ctj and (3j are arbitrary nonegative integers with a = a, + ■ ■ ■ + ag. 
This new Schwartz space is often denoted iS(M'') (observe that iS(M) = S). 
The properties of 5(]R'^) are discussed in standard textbooks [871 [90] give 
that the triplet 



where >S^(M'?) is the antidual of S{W), is a new RHS. 

The Schwartz space in all its versions has a topological advantage over 
other topological vector spaces: it is metrizable, i.e., its topology can be 
derived from a distance compatible with the structure of vector space. There 
are other interesting topological features on iS(]R'^) such as nuclearity that 
may be interesting like, for instance, the Dirac representation theory for 
some operators |3l|52], although this properties will not be discussed in here. 
One of the interests of the Schwartz space is that can be used to construct 
other types of RHS useful for the precise definition of the Gamow vectors for 
resonances. We shall discuss this construction in the next subsection. 

Other type of RHS includes nonmetrizable test vector spaces. The most 
popular example in this category is 

V C L2(M) c V , 

where V is the vector space of functions in S that vanish outside of a bounded 
interval. The space V is dense in L^(M) (and hence in S) with the Hilbert 
space topology r-^ and in S with its metrizable topology. Therefore, in order 
to distinguish V from S, we have to endow V with a topology different from 
the inherited from S. This topology is nonmetrizable but still keeps some 
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properties that make it tractable almost as if it were metrizablj^. 

This type of RHS can be generalized to higher dimensions in the con- 
figuration space. The vector space V{MP) is the subspace of the functions 
iS(M^) that vanish outside a compact set in MP, endowed with a nonmetriz- 
able topology of the strict inductive limit type [100] . We he the following 
RHS 

V{W) c L^W) c I?(RP)^ , 

for each positive integer p. Note that for p = 1, we obtain the previously 
discussed case. 



A. 2 Rigged Hilbert Spaces of Complex Analytic Func- 
tions. 

So far, we have introduced the most basic examples of rigged Hilbert spaces 
of functions. We shall not extend on their properties into here, addressing 
the interested reader to the speciahzed literature [SSI EZ]- With the help 
of those, we can construct the triplets relevant for the definition of Gamow 
vectors. 



A. 2.1 Rigged Hilbert Spaces of Hardy functions. 

To begin with, let us denote the open upper half plane (the set of complex 
numbers with positive imaginary part) as C"*". Analogously, C~ will be the 
open lower half plane, containing the complex numbers with negative imag- 
inary part. 

A Hardy function f+{z) on the upper half plane C"*" is a complex analytic 
function on C+ such that [Ml gSl EH [65] : 

/oo 
\ f {x + iy)\'^ dx < K < oo , (A. 11) 

-oo 



where K is certain positive real number. The inequality (lA.lip means that 



for each y > the function on the variable x, f{x + iy), is square integrable 
and all the integrals for y > are uniformly bounded by K. 



-•^^This is a strict inductive limit topology[100J of metrizable spaces. Problems of con- 
tinuity of operators on V can be investigated in each one of the metrizable spaces that 
construct the inductive limit, because of the Dieudonne-Schwartz theorems [100^ . That is 
why we said that the topology is "almost metrizable" . 



89 



The set of Hardy functions on the upper half plane form a vector space 
that we denote as "H^. Functions in "H^ have well defined (save for a set 
of zero Lebesgue measure) boundary values on the real line R. The func- 
tion of these boundary values is also is also square integrable. The 
Hardy function f+{z) determines uniquely (save for a set of zero measure) 
the boundary function /+(x). The reciprocal is also true and is a conse- 
quence of the well known Titchmarsh theorem [HT], so that /+(a;) determines 
uniquely the Hardy function f+{z) and, therefore, we can identify both. 

The definition of Hardy function on the lower half plane is similar, chang- 
ing y > hj y < 0. The vector spaces of Hardy functions on the upper and 
lower half planes are denoted by "H^ and "H^ respectively. Relevant proper- 
ties of these spaces can be found in the standard literature [Ml 1121 EH ES] • 
We list a few of them in here: 

i. ) It is very convenient to identify each function in f±{z) G with 
the function f±{x) of its boundary values on the real line R (allowed by 
the Titchmarsh theorem that also applies for "H^). The functions f±{x) are 
square integrable and therefore C L^(M). 

ii. ) There is a simple way to construct both Ti^ and T-t^ which is a 
straightforward consequence of the Paley-Wienner theorems [SOI ESI 1^ . 
which state that the Fourier transformation J-" is a unitary mapping between 
the Hilbert spaces L^(R'f) and "Hi, i.e., 

J^:L2(R+)=H^ ; : L^{Mr)=n\. (A.12) 

iii. ) Consequently, T-L^ are closed subspaces of L^(R). Furthermore: 

L\M.)=nl®Vl. (A.13) 

Thus, the Hilbert space of square integrable functions on the real line is a 
direct orthogonal sum of the spaces of Hardy functions "H^ and "H^. Then, 
any square integrable function f{x) can be uniquely written as 

/(x) = /+(a:) + /_(x) 

with /±(x) e n\. 

iv. ) Hardy functions can be recovered from its boundary values on a 
semiaxis of the real line. As a consequence of a result by van Winter [99] 
all values of any function in 1-L\ (including boundary values on R~) can be 
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obtained from its boundary values on the positive semiaxis M^, by means of 
a formula that uses the Mellin transform ^5] . 

V.) All the above properties (and some further results) are essential in 
the construction of new rigged Hilbert spaces [211 |5T] : 

nlns cnlc {nlnsy , (a.i4) 

where S is the Schwartz space. The space T-L^. H S is endowed with the 
topology on S and is a closed subspace of S. 

The van Winter formula establishes a one to one correspondence between 
Hardy functions and the functions of their boundary values on M"*". Thus, 



if we denote by "H^ H S 
have that 



the space of these boundary value functions, we 



nlns 



C 



nins 



(A. 15) 



is a pair of new RHS. The spaces Hi fl S 



are endowed with an special 



topology which is derived from the topology on S [SH [23]. 



A. 2. 2 RHS of entire analytic functions. 

Let us consider the vector space of the Schwartz functions that vanish outside 
a compact set included in the negative semiaxis (these functions are zero on 
the positive semiaxis M"*"), here denoted as ©(M^). The space of the Fourier 
transformations Z_ = J-'(P(]R~)) is a vector space of entire analytic functions 
which are Hardy on the upper half plane. Analogously = J-'(P(R+)) is a 
space of entire analytic functions which are Hardy on the lower half plane. 
The spaces Z^ have a nonmetrizable topology obtained from the topology 
on V{M.) [56n . We have new RHS given by 



and 



z^cnic zi 



(A.16) 



C [Z. 



(A.17) 



^"The spaces Z± are not discussed in this reference, although their properties are obvious 
from those of Z = J^(I?(R)) discussed in the given reference and the properties of Hardy 
functions. 
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These RHS have not been of great use in the description of resonances, 
although they are of some interest in the discussions of irreversibihty or time 
reversal invariance of the current formalisms on resonance phenomena |35] . 

A. 3 A rigged Hilbert space for the Friedrichs model. 

In the case of the standard Fridrichs model as presented in Section 2, we can 
easily construct the rigged Hilbert space in which all the generalized objects 
that appear in this situation make sense. 

The point of departure is the free Hamiltonian Hq in fl2.9p . According 
to the Gelfand-Maurin theorem, introduced in Appendix A, there is a rigged 
Hilbert space ^> C 7/ C such that i.) Hq<^ C ii.) Hq is continuous on 
$ with its own topology and iii.) for any uj G M"*^, the absolutely continuous 
spectrum of Hq, there exists a |ci;) G with Hq\uj) = uj\u). 

Let us recall (lA.Sp . let us omit the arbitrary Lp^ip & ^ and let us choose 
f{uj) = u) and A = Hq (in order to be consistent with the usual notation for 
the eigenvalues of Hq, we shall use cu instead of A in flA.7|) and (IA.8|1 ). Then, 



according to the Gelfand-Maurin theorem, one gets for the continuous part 
of Hq the following generalized expansion 

\u){uj\dfi{u) . (A.18) 







As the continuous spectrum of Hq is purely absolutely continuous and non- 
degenerated, (lA.lSP can be written as 



\uj){uj\du. (A.19) 

In addition to the continuous spectrum, Hq has an eigenvector that we have 
called |1) {Hq\1) = uq\1)). Then, the complete expansion of Hq in terms of 
its generalized eigenvectors with eigenvalues on its Hilbert space spectrum is 
given by f l2.9p . i.e., 

POO 

Hq = uq\1){1\+ \u){u\du. (A.20) 
Jo 

The explicit form of a RHS that implements in our case the conditions 
of the Gelfand-Maurin theorem is easy to find. In fact, we shall have a pair 
of RHS that satify these conditions. In addition, both members of the pair 
of Gamow vectors that describes the resonance will be functionals on each 
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of these RHS. These RHS are very easy to construct. In fact, as has a 
bound state and a non-degenerate continuous spectrum covering the positive 
semiaxis, the minimal Hilbert space in which Hq acts is given by 

n = C®L\R+), (A.21) 

where C denotes the set of complex numbers and © direct sum. All vectors 
in (OH) have the form fl^ . 

In our case, the vector space of test functions should have certain ana- 
lyticity conditions so that the antiduals contain the Gamow vectors. Several 
motivations (see [2llil|23l|29l|25l[271[26lliT]) suggest the use of functions 
which are at the same time Hardy and Schwartz, so that the spaces of test 
vectors are given by 



c^lsnnl ) , (A.22) 



where S denotes the Schwartz space, Ti^ the space of Hardy functions on 
the upper (with +) and lower (with — ) half planes. These functions are 
restricted to the positive semiaxis because they represent wave functions in 
the energy representation and we assume that the energy is always positive. 
Note that the values of a Hardy function on the positive semaixis determine 
all its values as a consequence of a Theorem due to van Winter [99]. Finally, 
the RHS for the basic Friedrichs model is the following: 



(A.23) 



The 



Note that the duals $± = 

Gamow vectors n = 0, 1, belong to and the Gamow vectors \ fn) 
n = 0, 1, belong to 

We recall that the potential is given by ( I2.10p . i.e.. 



V 



{r{co)\co){l\ + fiu)\l){u\)doj. 



(A.24) 



Note that if f{co) is a square integrable function, then y is a bounded 
operator from % into "H. The proof goes as follows: Take ip as in (12. 2p . then 
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< 



2 poo 
^0 



'iiy.ir+i«r"^"' 



(A.25) 



Since is arbitrary, we conclude that V is bounded. In addition, V is 
Hermitian. Then, a consequence of the Kato-Relhch theorem [87] estabhshes 
that H := Hq + V self adjoint on "H. This is true for any square integrable 

form factor /(w). If, in addition, /(w) G 5 fl 1-L\ , H is continuous on $| 

respectively and can be extended into a weakly continuous operator on the 
antiduals. Then, the action of H on the Gamow vector makes sense. 

It could happen that f{uj) is not in this class, still being square integrable, 
or even that /(cj) were not square integrable, being for instance a polynomial 
in the variable u. In this cases and in some others (like for instance if the 
form factor /(cu) is a Dirac delta), if is a continuous mapping from the space 
of test functions into the antiduals: 



H : $± 



(A.26) 



In some cases, the domain of H can be extended to some subspaces of 
like for instance to the spaces spanned by the Gamow vectors. This happens 
for instance in the two situations mentioned before of square integrable or 
polynomial form factor, also the form factor ^/uj/P{u) for the double pole 
resonances in the Friedrichs model is of this form. Then, the extension of 
H into the antiduals cannot be produced by the duality formula ( lA.Sp . but 
instead by completion [16j. Note that this type of extensions may not be 
extended into the whole antiduals. 

Other types of more sophisticated Friedrichs model need different RHS 
construction. Along the present work, we make a case by case study of these 
different types of constructions. 
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B Calculation of 77(2:) in the ordinary Friedrichs 
model. 



In this appendix, we want to obtain the expression of the reduced resolvent 
for the simplest form of the Friedrichs model leading to equation (12.141) . We 
shall make use of operator theory and follow the guidelines of jSj. In fact, 
we want to show that 

^ -11)= f-z + Uo + X' r^-^duY'. (B.l) 



z — H \ Jo z — oj 

We start with the demonstration of the so called Friedrichs condition that 
reads 

QVQ = 0, (B.2) 

where O is the zero operator. The action of O on any vector ip & is the 
zero vector G "H. Taking ip ETi as in (12. 2p . we have that 



QVQ^l! = QVQ 



a 




g = . (B.3) 








Thus, flB.2p is proven. 

The second step in the proof of fIB.ip is noting that the projectors P and 
Q commute with the free Hamiltonian Hq. Since Q = I — P, this claim 
would be shown if we can see that PHq = HqP. In fact, if the arbitrary 
vector e ?^ is as in (12. 2p . we have 



pho^ = PHo 



a 






a 

HoPtjj = HoP { ]=Ho\ 1 = 1 1 , (B.4) 
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as we claimed. 

Now, we use the second resolvent identity as given in |101] : 

R{z, H) = R{z, Ho) - X R{z, Ho) V R{z, H) , (B.5) 

where 



R{z, H) = = {zl - H)-' , R{z, Ho) = ^ (zl - Ho)-' 

z — H z — Ho 

and H = Ho + V . This gives: 

P^-P = P^—P-XP^—V^-P. (B.6) 
z — H z — Ho z — Ho z — H 

Now, if / is, the identity operator, as usual, we can insert P + Q = I m the 
last term of (]B.6P as follows: 



+P—^PVQ-^P + P—^QVQ-^P. (B.7) 

z — Ho z — H z — Ho z — H 

As P and Q commute with Hq and PQ = O, we have that 

P^—Q = 

z- Ho 

and therefore the second and forth terms in the right hand side of (IB. 70 
vanish. We also have that 

z — H z — Ho z — Ho z — H 



The first term of the right hand side of flB.7p also vanishes. Then, we have: 
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z — H z — Hq - II 

Z — Hq Z — H Z — Hq 

-\Q—^PVQ-^P-\Q—^QVP^-P (B.9) 

Z — Hq Z — H Z — Hn Z — H 



z — 


H 


1 




z — 


H 


1 





Again, the second and the third term in the right hand side of ( IB.OP vanish. 



Due to the Friedrichs condition (]B.2p . also the first term vanishes. Thus, 



Q^—P=-\Q^—QVP^—P (B.IO) 
z — H z — Hq z — H 

Inserting (12.91) into (12. 6p and then, (12. 6p into (12.51) . we have 

P^-^P = P^^P-XP^^PVP^-P 

Z — H Z — Hq Z — Hq Z — H 

+X'P^^VQ^—QVP^—P. (B.ll) 

Z — Hq Z — Hq Z — H 



If we multiply fIB.lip to the left by {z ~ Hq) P, we can write 

iz-HQ)P-^P = P 

Z — Uq 

and therefore, 

P = [{z - Hq) P + X PV P - X'' PV Q QV P] P ^— P , (B.12) 

z — Hq z — H 



or 



where 



P P = [G{z)]-' P , (B.13) 



G{z) = [{z - Hq) P + XPV P - X^ PV Q — ^-—QVP]. (B.14) 

Z — Hq 
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Since, 



PVQ^—QVpI ^ I =PVQ^—Qv( ^ 







PVQ 



Z-Hr 



Q 



PVQ 



z-H> 



V /*M 



PVQ 



\ / I/Ml' du 



P 



2 — W 



(B.15) 



and P y P |1) = 0, we have that 



G{z)\l) 



Consequently, 



z — (jjq — y 



z — u 



duo 



(B.16) 



z-H 



z-H 



which proves fIB.ip . 



,2 r i/MP . ' 

Z — Wo — / 



(B.17) 
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